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Abstract 


-^•Simple  models  of  polymer  dynamics  are  available  in  dilute  solution, 
moderate  concentrations  and  melts,  since  it  is  possible  to  make  models 
of  the  motion  in  these  cases.  A  series  of  power  laws  result  which  fit  well 
with  computer  simulation.  It  is  more  difficult  to  derive  these  models 
directly  from  sensible  equations  of  motion,  but  progress  in  this  direction 
is  reported  in  the  paper.  - 
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I  §  |  Introduction 

Polymer  dynamics  involve  different  time  scales  according  to  the  part  of 
the  polymer  involved  and  the  environment  that  it  finds  itself  in.  In  smaller 
j  molecule  dynamics  one  can  identify  three  obvious  regimes,  the  molecule  alone, 

usually  meaning  in  the  gaseous  phase  where  one  has  an  n  body  problem  if 
there  are  n  atoms,  an  n  body  problem  with  external  noise  and  friction  which, 

(  roughly  speaking,  is  what  one  expects  in  a  liquid,  and  to  some  extent  a  glass, 

and  finally  part  of  an  Nn  body  problem  in  a  crystal.  Notoriously  the  middle 
problem  is  the  worse  since  the  first  has  non  stochastic  equations,  and  the 
I  crystal  is  at  least  quite  explicit  >  and  there  is  a  well  defined  band  theory. 

The  liquid  case  is  the  most  difficult  since  the  molecule  has  a  structure 
comparable  to  its  surroundings,  it  could  for  example  be  in  the  liquid  phase  of 
j  its  own  species,  and  there  is  no  easy  way  to  simplify  the  environment  as  there 

is  with  a  large  Brownian  particle  in  a  small  molecule  fluid. 

If  one  studies  a  polymer  which  is  a  string  of  smaller  monomers  then  only 

*  those  aspects  of  its  behaviour  which  have  to  do  with  the  polymerized  nature 
have  any  hope  of  being  easier  to  understand  than  the  corresponding  monomeric 
liquid,  and  the  solid  state  of  polymers  will  always  be  more  difficult  than  the 

*  crystalline  phase  of  the  monomer,  although  again  they  may  have  some  features 
which  are  simpler  than  the  monomeric  glass. 

One  can  consider  the  behaviour  of  polymers  characterized  by  relaxation 
times  or  typical  frequencies,  and  the  discussion  above  amounts  to  saying  that 
when  one  studies  high  frequencies,  a  polymer  melt  is  only  marginally  different 
from  a  monomeric  melt,  and  indeed  this  is  borne  out  experimentally.  But  the 
polymerized  nature  of  the  material  totally  alters  melt  and  solution  behaviour  at 
low  frequencies,  the  long  chain  nature  being  reflected  in  the  spectacular 
increases  in  viscosity  which  are  so  characteristic  of  polymers.  There  will  of 


i 
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course  be  effects  in  the  solid  state  also,  but  they  are  not  so  spectacular 
since  the  fall  of  temperature  means  that  Van  der  Waals  or  hydrogen  bonds  are 
now  real  bonds  and  all  glasses  are  polymer  networks  in  some  sense;  a  glass 
made  of  polymer  is  much  the  same  as  any  other  glass  once  one  is  well  past 
the  glass  temperature. 

So  if  one  wants  to  study  stimulating  theoretical  ideas  it  is  natural  to 

aim  at  effects  which  are  characterized  by  long  time  scales,  and  an  obvious 

pitfall  ensues.  To  graft  that  kind  of  theory  out* monomeric  structure  is  a 

hard  problem,  but  one  inescapable  if  that  is  what  one  really  wants  to  do; 

it  is  no  good  solving  by  methods  which  are  fine  on  long  time  scales  and 

extrapolating  to  short  time  scales.  This  is  a  general  caution  to  be  noted. 

Current  polymer  dynamics  initiated  in  the  famous  Rouse  paper  is  at  the  level 

of  elastic  wave  theory  of  solids,  the  Bom  Von-Karman  theory,  or  equivalently 

it  is  equivalent  to  the  Rayleigh,  Jeans  and  Wien  formulae  of  black  body  theory. 

The  Bom  Von-Karman  theory  does  not  produce  that  multiplicity  of  branches  in 

the  vibrational  spectrum,  and  in  the  way  that  faults  affect  it.  But  there 

is  an  enormous  richness  in  polymer  theory  at  the  Rouse  level  since  a  whole 

host  of  new  problems,  still  entirely  confined  to  long  time  scales,  appear  in 

'  polymer  problems,  and  they  pose  problems  which  are  very  deep  mathematically 

but  intuitively  obvious ,for  example:  how  long  will  it  take  for  a  knot  to  untie 

if  the  string  is  subjected  to  Brownian  motion?  (If  the  string  is  in  a  melt, 

/  the  melt  encloses  it  like  a  tube,  and  one  question  is  asking  how  long  it  takes 

for  the  chain  to  wriggle  out  of  this  knotted  configuration  into  a  new  configuration. 

This  is  clearly  a  shorter  time  than  wriggling  right  out  of  the  initial  tube  into 

■i  an  entirely  new  tube,  and  that  time  can  be  calculated  by  reptation  theory  and 

-2 

is  proportional  to  L  .  So  unless  the  new  tube  is  the  same  class  of  knot  one  has 
the  solution  to  the  query  at  a  crude  level). 
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I  shall  put  forward  a  chain  of  developments  of  dynamics  in  recent  years 
which  amount  to  a  reasonable  intuitive  solution  for  long  time  scales,  and 
point  out  the  basic  problem  of  improving  these  theories  to  the  level  of 
mathematical  proof. 

§2.  The  single  chain 

Forces  between  monomers,  and  also  forces  between  the  polymer  and  its 
surroundings  have  the  effect  of  expanding  or  contracting  the  chains.  This 
paper  is  concerned  with  dynamics,  so  it  will  omit  all  discussions  of  swelling 
or  precipitation.  This  is  a  topic  'under  control'  and  can  be  built  into  the 
analysis  without  too  much  difficulty;  it  is  not  worth  while  muddying  the 
equation  with  complexity  which  will  not  be  used.  There  is  a  dynamical 
equivalent  to  this  simplification  in  that  every  motion  involves  crossing 
potential  barriers  in  quite  complex  ways,  and  these  ways  are  affected  by  the 
environment  just  as  a  single  particle  can  have  complex  jump  behaviour 
but  still  in  the  end  will  have  its  long  time  motion  governed  by  a  simple 
diffusion  equation; 


(2.1) 


one  can  hope  for  a  similar  equation  for  the  polymer. 

In  the  presence  of  a  potential  V(x)  a  particle  in  equilibrium  will  have 

P  .  < 

o 

where  is  the  normalisation.  This  modifies  (2.1)  to 
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a  simple  but  rather  unfamiliar  form  since  one  generally  does  not  find 
Brownian  particles  in  central  fields  of  force.  The  mathematical  form  is 
familiar  in  the  Fokker-Planck  equation  for  the  case  of  particles  with 
inertia  and  soft  forces  which  has  the  form: 


.7}  ^V)being  the  velocity  dependent 
tensor  diffusivity  required  in  say  plasma  physics. 


o 

(2.3) 


What  is  Pq  for  a  polymer:  In  equilibrium  it  is  a  random  walk.  Suppose 

it  were  freely  hinged  rods  whose  ends  were  at  Yt  ,  _  and  each  rod 

is  of  length  : 


■'J  S'lT;  -Tu,  /  -i) 
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(2.4) 


It  is  well  known  that  the  end  to  end  statistics  or  indeed  the  statistics  of 
any  points  distant  from  one  another  along  the  chain  is  independent  of  the 
precise  neighbour  configuration,  and  the  analogue  of  the  variable  w&L 
describ&selastic  waves  in  a  solid,  or  the  macroscopic  density  or  velocity  in 
a  liquid,  are  the  fourier  components 

R r 

(2.5) 

which  gives(by  a  well  known  straightforward  derivation  not  reproduced  here) 
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(I  have  used  the  complex  form  for  simplicity,  but  it  contains  the  real 
form  if  required  by  relations  between  the  real  and  imagining  parts  of  R^  ) 
The  q  in  the  sum  cannot  exceed  the  number  of  links  (recall  the  Born  theory) 
but  in  fact  if  ever  we  need  to  restrict  the  sum  the  theory  will  be  wrong 
since  high  q  brings  us  into  a  region  outside  the  present  validity.  Note  that 
there  is  no  kT  in  P  ,  the  random  walk  has  entropy  but  not  fins  internal 
energy.  The  analogue  of  (2.2)  is  now 
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Provided  one  realises  that  only  the  q  variables  have  meaning^one  can  also 


back  Fourier  transform 


-  v/V  e 


^  rL^)^s 


2  _  Us  2  p  '2  -  3&y  ' 

9t  j  ^5^)  i  y. 


p= 


(2.8) 


where  the  continuous  arc  lengthjs  have  replaced  the  discrete  ,  but  for  those 


who  find  functional  differential  equations  hard  to  swallow  one  may  emphasize 
that  (2.7)  is  exactly  the  same  as  a  equation  for  highly  damped  phonons  in  a 
solid,  or  radiation  in  a  highly  emitting  and  absorbing  material. 

Equation  (2.7)  is  the  Rouse  equation  and  the  q’s  label  Rouse  inodes.  It 
can  be  solved  exactly,  for  if  one  writes 


p  =  p  *'/&,!*) 


one  has 


y  /_?1-  -t  P  ,*/«,/*  ~  ^ 

2.  i  ^  ^  ,A  If  ■»%'  I+Hi  ~ 


f  \  1 5ft  3ft  , 
M  ~  7 


,  \\P=o 

al-Srt.  <2-10> 
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which  is  Hermites  equation  and  is  in  all  the  books  under  the  quantum  theory  of 
the  harmonic  oscillator,  give  or  take  a  few  constants  and  i. 

From  this  equation  one  can  work  out  the  probability  of  the  polymer  getting 
from  any  given  shape  to  any  other.  By  taking  the  first  moment  of  (2.7)  we 
see  that  it  is  derivable  from  a  Lange  min  equation; 


-v  R 


3 <T 

cJt xl 


= 


(2.11) 


where 


is  a  friction 
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(2.12) 
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'•  V 

0 

V* 


is  the  free  energy  per  mode, 

X  being  a  random  force.  This  friction  V  R,  is  just  like  the 
f  q  a  Brownian  sphere  in  Stokes  formula,  which  amounts  to  sayi 

that  the  chain  suffers  a  Stokes  drag  along  its  length  and  this  drag  shows  no 


coupling-act*^  c.Ka.0^  . 

In  fact  it  must  show  coupling  since  the  flow  around  any  point  will  affect 
neighbouring  parts  of  the  chain,  unlike  Stokes'  problem  where  each  Brownian 
sphere  is  considered  remote  from  each  other  one.  The  hydrodynamic  effects  have 
been  put  in  by  Ziroro  and  modify  the  diffusion  equation  in  a  way  most  easily 


written  in  the  form  (2.8). 

r  ?  rfr: 
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where  e/.  f ^  label  a  pair  of  polymers,  ^  are  Cartesian  indices 

and  2)  is  known  in  the  hydrodynamical  literature  as  the  Oseen 

n 

tensor.  It  is  most  easily  written  in  fourier  transform 

i'*-)  -  V"  -  ^  K/k 1 

(The  speed  of  propogation  of  the  signal  through  the  liquid  is  here  ignored). 

Equation  (2.10)  is  frequently  simplified  by  replacing  D  by  its  average  for 
a  random  walk: 
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(2.15) 
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so  that 


becomes 
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From  these  forms,  or  directly  from  the  Langevin  equation  one  can  derive 


the  correlation  functions 

0))y 

=  e-  (2.ia) 

.  s^yVj  (,w 

for  the  Rouse  case  and  ^  ^ 

f  a,  0  -  '>s)  /^y'.  /2 

,  p*j£)*  i  i  J 

(2.20) 

for  the  Zimm  case.  The  Zimm  result  is  more  realistic  since  it  is  derived 
from  possibly  realistic  model,  that  of  a  polymer  embedded  in  a  fluid  with 
normal  hydrodynamics.  The  weakness  of  the  model  lies  in  the  neglect  of  the 
entanglements.  The  Rouse  model  on  the  other  hand  is  not  realistic  as  it 
stands  since  there  is  no  justification  for  assuming  a  drag  coefficient  in  the 
form  of  a  simple  constant,  which  certainly  will  not  result  from  the  equations 
of  hydrodynamics.  Experimentally  the  viscosity  of  a  dilute  polymer  solution 
verifies  the  Zimm  result,  for  the  viscosity  in  this  case  is  proportional  to  [_ 

(i.e.  where  M  is  the  molecular  weight). 
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screened  once  the  density  of  polymer  is  such  that  the  chains  overlap.  If 
one  takes  equation  (2.14)  and  averages  away  all  polymers  except  one,  one 
finds  that  simple  screening  does  obtain,  and  just  using  mean  densities  the 
equation  governing  the  effective  hydrodynamics  equation. 


du 

d-t 


t  v  fr)  -f  Yl 


and  the  effective  polymer  dynamics 

+  *7  f  ^  7rf) 


(2.21) 
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V  (2.23) 
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If  an  arbitrary  friction  term  is  still  left  in,  it  will  appear  as 


•J  -4-  v  Jcp  )gy  -h  =  21 


(2.23) 
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At  high  enough  densities  one  can  expect  a  solution 


(Vjg  4-  $  \j 


where,  by  studying  '  2  •  finds  that 
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=  y}*(v-+ $~y  l‘ili  >/ 


(2.26) 


(2.27) 


This  is  of  course  an  asymptotic  solution  to  have  so  simple  a  form 
and  the  full  equations  are  quite  complicated.  A  particular  problem  is 
that  the  entanglements  must  dominate  at  high  enough  density,  and  it  is  not 
clear  if  this  screened  regime  is  ever  realistic;  but  experimentally  one 

certainly  has  a  regime  of  viscosity  M  before  the  melt  behaviour  characterized 

k  »  3  + 

by  /''j  sets  in. 

This  region  will  be  considered  below. 

The  conclusion  of  this  section  is  then  that  the  diffusion  of  a  single 
point  on  a  polymer  in  solution,  or  the  two  body  correlation  function,  can  be 
obtained  provided  that  entanglements  are  ignored,  from  the  formula 


?  f~  I 


'  —  e  (o^28) 


(one  point  function  has  S  «  S  ,  2  two  body  correlation  integrates  over 

s-s') 

where  J  is  given  by  (2.24)  coupled  with  (2.23). 
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£ 3.  Entanglements 

We  have  seen  in  the  last  section  that  the  dynamics  of  chains  which  are 
transparent  to  themselves  can  be  resolved  which  they  are  in  Brownian  motion 
with  a  simple  drag  friction,  or  surrounded  by  a  viscous  fluid,  and  when  many 
chains  are  present  some  screening  phenomenon  intervenes.  But  in  fact  all 
chains  are  non  transparent,  and  this  fact  dominates  motion  in  melts  or 
concentrated  solutions  since  almost  all  motions  are  blocked.  A  surviving 
motion  will  be  the  wriggling  of  the  chain  up  and  down  the  statistical  tube 
which  one  can  imagine  surrounding  it  on  average,  and  a  plausible  visco  elastic 
theory  can  be  developed  from  the  reptative  motion.  It  is  also  possible  that 
some  kind  of  collective  motion  could  exist,  but  no  convincing  theory  of  this 
exists  at  present.  To  illustrate  the  concept  of  reptation  there  is  a  simple 
computer  experiment. 

Suppose  a  lattice  is  made  up  in  two  or  three  dimensions  which  we 
illustrate  by  a  two  dimensional  array  of  dots.  The  polymer  is  moved  by 
Monte  Carlo  moves.  It  is  well  known  that  without  the  lattice  this  gives  the 


Rouse  equation.  If  one  regards  the  lattice  as  defining  a  tube  one  can 
characterize  this  by  pulling  the  polymer  taut  one  generates  another  walk 
with  a  larger  steplength,  called  the  primitive  path.  The  steplength  of  this 
primitive  path,  a  say,  is  also  the  thickness  of  the  tube. 


:  y-v 
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Suppose  Che  distance  along  Che  Cube  is  S  and  the  absolute  coordinate  R 

iV  =  a[S,  ~S^j 

(3.1) 


Now  consider  a  point  on  the  polymer  initially  labelled  and  after  a  time 


labelled  .  Then 

y  _ _  S  ^  ^ 

>  *■  '  '  is  given  by  the  Rouse  equation  in  one  dimension 

and  this  as  was  shown  above  gives 
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V 

J  t- 

O) 


C&  Q  S' Oka. 


(3.2) 


(Note  the  distinction  between  this  and  a  Brownian  point  particle  which 
gives  b  ) . 

<r  fefir)  -  S(o))  >  OC  Jt 

y  (3.3) 


R(s,o)^y 

PC  £,lit‘ 


(3.4) 


A 


Thus  one  might  expect  a  time  scale  initially  of  t2  for  the  polymer  has  yet 
to  notice  the  tube.  Then  t^,  it  is  diffusing  up  and  down  the  tube.  The 


centre  of  mass  of  the  time  is  the  special  component 

.1 


?  = 


=  O  and  will  have 


the  usual  t  for  S  and  hence  t2  for  R,  so  after  a  longer  time  again  the  c.m. 


will  diffuse  and  take  any  point  with  it,  gfing 


Finally  the  polymer  gets  right  out  of  the  tube  and  creates  a  new  tube 
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and  the  final  t  lav  must  result.  Thus  one  expects  the  hierarchy 


(3.5) 


The  initial  t^  is  a  bit  difficult  to  find  with  a  computer  which  is  moving 
the  system  very  quickly  indeed  to  get  into  the  other  ranges,  but  with  the 
eye  of  faith  one  can  see  the  other  three  ranges. 


9 


t 


9 


9 


* 


f 


2.0  3.0  4.0  S.O  60 


This  is  for  a  fixed  lattice  and  a  free  polymer.  When  all  the  polymers  move 
the  time  taken  over  the  computation  has  to  be  immense  both  because  one  has 
to  move  many  more  polymers  and  the  checking  of  the  topological  integrity  is 
much  more  difficult  and  no  longer  can  be  incorporated  into  the  program  without 
the  use  of  memory  and  because  the  system  becomes  'soggy'  and  collective  modes 
are  clearly  possible,  but  not  easy  to  incorporate  in  a  one  point  Monte  Carlo 
system.  Various  authors  give  various  results,  but  my  belief  is  that  the 
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reptation  picture  has  to  be  the  starting  point  of  any  theory.  The  authors 
papers  with  K  E  Evans  illustrate  a  whole  range  of  dynamical  effects  of  the 
fixed  lattice  models,  and  these  results  have  been  extended  by  R  J  Needs  to 
star  molecules,  verifying  the  results  of  the  theories  of  de  Gennes  and  of 
Doi,  and  the  experimental  diffusion  results  of  J  Klein. 

Thus  we  can  now  argue  that  the  diffusion  law  as  a  function  of  concentration 
will  go  like 

(hydrodynamics)  2/3 

Low  concentration  (  Zimm  )  t  finally  t 


Higher  concentration  (Screened) 

(to  Rouse) 


t2  finally  t 


1  1  1 

(Reptation:  de  Gennes)  t  ,  t*,  finally  t2  then  t. 


The  intermediate  regions  are  difficult  and  a  principal  difficulty  is 
that  the  fundamental  problem  of  how  to  describe  the  motion  of  entanglements 
has  to  be  resolved  in  the  region  between  higher  concentration  and  very  high 
concentration.  In  the  former  (and  in  the  Zimm  regime)  one  just  ignores  it. 

In  the  very  high  density  regime  it  becomes  straightforward  once  it  can  be 
represented  by  a  tube. 

How  can  one  give  a  representation  of  the  tube  in  mathematics?  One  way 
is  to  consider  the  polymer  at  two  times.  Since  this  paper  is  being  read  at 
an  NMR  conference  I  will  not  worry  about  very  long  time  scales;  alternatively 
we  can  consider  rubbers  where  the  tubes  are  permanent.  Suppose  the  polymer 
is  initially  R2(S)  or  R2(q)  and  later  at  R, ((J) .  Then  over  the  time  scale  over 
which  the  polymer  moves  a  tube  step  length  it  will  not  encounter  the  walls. 
(This  whole  argument  is  going  to  be  ’average’).  Let  this  time  be  called^ 


where  <^h as  the  dimension  of  q  and  is  related  to  a,  'j 0  ^  to. 


(3.7) 
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The  Rouse  equation  is  a  version  of  Hermite’s  equation  and  the  joint 
probability  of  finding  at  t=o  and  R2^)  at  t  =  t  is  given  by 

)=7T i1  Aa*)1-) 

•  (foil1*  18^ I*)  eHUIXj’t  +  <JJ  /P,,  Bj,  +5i-) 
xf(WA.  ]>  1^0 /»*>•+■  j’3'1  <3’8) 


Note  that 


p  %.7,r6.-],t)7r  <*«, 

; 


where 


P  -PurV 

P  -  (/<*„  A  3  P,  P  ^ 

1  »  -  t  il  arr  t  ’  y . 


(3.10) 


This  is  for  a  free  polymer. 

We  now  argue  that  for  a  polymer  in  a  pipe  the  joint  probability  is  the 
expression  above  with  "£) 

"  PCEjgd,  fej  ^ 

Note  the  complexity  of  the  joint  form.  It  would  have  been  quite  wrong  to  solve 
the  problem  of  a  single  polymer  constrained  to  one  dimension.  The  dynamics 
of  a  polymer  R  tied  by  a  tube  constraint  represented  by  a  locus  R2  can  be 
developed  by  noting  that  is  the 

equilibrium  distribution  of  R  hence  the  equation; 


3 

at 


r>  ^  (} 


M 


*■ 


It  will  be  seen  from  the  structure  of  PQO  that  R  and  R£  are  strongly 
correlated  over  J  ^  0  if  S  ^  Q 


but  weakly  correlated  for 


7  >  <?«  <>.  s<  a.  ( 


which  is  precisely  the  primitive  path  picture. 
It  will  be  seen  that  — - —  9  Poo 


R 


tt  CK 


i .fi-Av"  -lrKpV4?Sit^ 


R  and  R^,  so  that  (3.11)  is  still  a  version  of  Hermites  equation  and  can 
still  be  written  down  in  closed  solution  so  that  one  can  now  derive  a  form 


for  the  probability  that  a  chain  constrained  by  a  tube  to  R2  starts  at 


and  is  at  R  at  t. 

But  the  answer  is  sufficiently  algebraically  tedious  that  I  will  not 
write  it  down  (it  possesses  the  property  of  being  tedious  even  though  trivial!) 

A  phenomenology  is  now  complete,  but  what  of  the  mathematics  of  entanglement 


4.  Arigorous  theory 

This  can  be  from  the  full  hydrodynamic  form  (  £ • ’H~ 

because  this  equation  comes  from  the  Stokes  boundary  condition 

£  -  &(£) 

where  is  the  fluid  velocity.  Hence  if  R.  »  R, 


I 
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It  is  now  possible  to  argue  that  this  condition  is  sufficient  to 
stop  the  polymers  crossing.  Hence  the  entanglements  are  already  in  (  2.' ^  )• 
A  particularly  simple  version  of  this  is  to  freeze  all  the  polymers  except  one. 
Then  one  can  show  that  the  last  polymer  ,  satisfies 


l 


So  at  a  rigorous  level  we  have  (4.1)  or  (  gtgjf})  which  'Vv'A *-*("  be 

equivalent  to  (3.11)  which  expresses  the  primitive  path  structure  explicitly, 
.  .  (2)  (3) 

i.e.  if  one  could  average  away  R  R  . it  would  have  to  be 

in  terms  of  one  of  the  primitive  paths  defined  by  these  loci.  The  point  is 
that  the  solutions  of  (4.1)  (as  shown  in  ref(  )  )  break  up  into  a  set  of 
probabilities  each  associated  with  one  of  the  primitive  paths  of  the  network. 


This  is  all  quite  straightforward  for  melts  for  in  spite  of  the 
apparent  complexity  it  is  merely  showing  a  pathway  from  basic  equations  to 
the  results  of  ^  3,  and  to  a  much  more  accurate  form  of  correlation  functions. 
I  have  hope  that  having  got  these  basic  equations  one  can  explore  the 
transition  region,  but  have  a  horrible  suspicion  that  computer  simulation  will 
get  there  first. 
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ABSTRACT 

A  procedure  of  analysis  of  NMR  measurements  is  proposed  as  an  approach 
to  the  observation  of  a  single  chain  diffusion  process  in  a  melt.  The  dynamic 
screening  effect  is  applied  to  the  definition  of  a  temporary  submolecule  ;  this 
is  used  as  a  semi-local  probe  to  investigate  collective  motions  of  all  parts 
of  a  polymer  chain.  Measurements  performed  on  "^C  nuclei  or  protons  lead  to 
similar  conclusions.  The  best  agreement  with  experimental  results  is  obtained 
by  combining  a  M3  dependence  of  the  terminal  relaxation  time  CM  is  the  chain 
molecular  weight]  with  a  multiple-mode  relaxation  spectrum  ;  all  modes  have 
the  same  statistical  weight. 
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NMR  AND  MACROMOLECULAR  MIGRATION  IN  A  MELT  OR  IN  CONCENTRATED  SOLUTIONS 


J.P.  COHEN  ADDAD 
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University  Scientifique  et  Medicale  de  Grenoble 
B.P.  68  -  38402  SAINT  MARTIN  O’HERES  CEDEX  (France) 


I  -  JHTR0VUCT10N 

""The  purpose  of  this  paper  is  to  analyse  the  migration  process  of  long 
polymer  molecules  in  a  melt  or  in  concentrated  solutions  as  it  may  be  observed 
from  the  dynamics  of  the  transverse  magnetization  of  nuclear  spins  linked  to 
these  chains. 

XThe  low  frequency  viscoelastic  relaxation  of  polymer  systems  is  known 
to  be  mainly  controlled  by  the  mechanism  of  dissociation  of  topological  cons- 
traints  excited  on  chains  and  which  are  called  entanglements ,  *  .  This  mecha¬ 
nism  exttbits  a  strong  dependence  upon  the  chain  molecular  weight. 


-These  topological  constraints  also  govern  the  diffusion  process  of  polymer 
chains  .  So,  the  accurate  description  of  the  diffusion  motion  of  a  chain  may 
be  a  convenient  way  to  characterize  disentanglement  processes  necessaffly  invol¬ 
ved  in  any  model  proposed  to  explain  viscoelastic  effects.  ■ _ 


Therefore,  it  is  worth  trying  to  define  an  experimental  procedure  giving 
a  direct  observation  of  the  diffusion  of  a  single  chain  at  a  molecular  scale 
although  it  is  in  dynamical  interactions  with  all  surrounding  chains. 


2.- 


Two  problems  are  encountered  In  attempting  to  observe  the  diffusion 
of  a  chain  at  a  molecular  scale. 

Space.- & cate,  of,  meaiuAementA 

One  of  these  is  a  problem  of  space-scale  of  measurement.  It  is  due  to  a 
dynamic  screening  effect  induced  by  topological  constraints  and  characterized 

O  ~ 

by  a  correlation  length  ag  ^  20  A.  Within  a  space  domain  defined  by  o^.  dyna¬ 
mical  fluctuations  concern  short  chain  segments,  only  (the  number  of  skeletal 
bonds  Ne  is  about  2  x  102)  j  they  correspond  t<*  high  relaxation  frequencies 
(-  10  Hz)  .  These  motions  do  not  depend  upon  the  molecular  weight  of  polymer 
chains.  There  are  no  correlations  of  fluctuations  of  these  short  segments  from 
one  screening  domain  to  another  one.  Therefore,  to  observe  a  significant  dis¬ 
placement  of  a  long  chain,  it  is  necessary  to  have  a  space-scale  of  measure¬ 
ment  longer  than  ae  j  otherwise,  only  local  motions  will  be  perceived.  This 
difficulty  could  be  overcome  by  performing  light  scattering  or  neutron  scat¬ 
tering  experiments.  But  there  is  also  a  problem  about  the  time-scale  of  mea¬ 
surement. 

The,  time.  bcaZe.  oj  meoAitnement 

The  time  interval  Tp  requires  to  observe  the  diffusion  of  a  chain 
over  its  own  dimension,  R^,  is  about  1  sec  for  a  chain  molecular  weight  M  v  io6 
and  Rg  ^  2x102  A  j  this  corresponds  to  a  diffusion  coefficient  Dseif  ^  10~10  cm2 
sec-1.  The  time  interval  is  about  equal  to  .1  msec  for  a  displacement  over  a 

O 

distance  smaller  than  Oe  =  20  A.  These  time  scales  are  not  appropriate  to  light 
scattering  or  neutron  scattering  experiments. 

Vlifax Aion  coeAAicicrvt  me/uuAmenti 

The  difficulties  about  time  scales  and  space  scales  have  been  overcome 
in  two  ways.  The  first  one  was  to  observe  the  diffusion  of  a  chain  at  a  macros¬ 
copic  scale  and  not  at  a  molecular  one.  Early  measurements  were  made  using 

chains  labeled  with  radio-tracers  Deuterated  chains  moving  through  a  pro- 

t  R 

tonated  matrix  were  observed  from  inrared  spectroscopy  0  ;  the  diffusion  coef¬ 
ficient  of  polyethylene  chains  measured  according  to  this  experimental  proce¬ 
dure  was  found  to  vary  as  the  square  of  the  chain  molecular  weight  within  a 
reasonable  accuracy. 


3.- 


Another  way  was  to  shorten  the  time  scale  of  chain  diffusion  by  consi¬ 
dering  short  chains  in  concentrated  solutions  or  long  chains  in  semi-dilute 
solutions.  The  diffusion  coefficient  of  polystyrene  chains  labeled  with  a  photo- 

chromic  probe  and  moving  through  a  pulsed  pattern  of  interference  fringes  was 

3 

shown  to  well  obey  the  predicted  formula  : 


0=  M-2^1*75 


C^(s  the  polymer  concentration  7 .  A  pulsed  field  gradient  NMR  technique  0  was 
recently  used  to  measure  the  diffusion  coefficient  in  the  concentration  range 
0.04  <  C  <  0.16  (W/W]  for  chain  molecular  weights  lower  than  3  x  105. 


It  will  be  shown  from  this  paper  that  quantum  coherence  properties  of 
nuclear  spins  may  be  used  to  investigate  a  chain  diffusion  process  not  by  mea¬ 
suring  any  diffusion  coefficient  D  but  by  adjusting  and  comparing  the  chain 
relaxation  spectrum  characterized  by  : 


*  0..,  V<R%> 


to  an  internal  NMR  reference  frequency  much  lower  than  the  Larmor  frequency. 


The  adiabatic  relaxation  of  the  transverse  nuclear  magnetization  M^Ct) 


nxct)  -  <Mx(o)  Mx C t ) > 


is  known  not  to  depend  upon  any  exchange  of  energy  between  the  spin-system 
and  the  thermal  bath  j  therefore  there  is  no  resonance  condition  to  be  fulfil¬ 
led  with  the  Larmor  frequency  ?  it  is  the  reason  why  slow  diffusional  proces¬ 
ses  can  be  observed. 


II  -  WO  BASIC  m  RESULTS 


Two  basic  NMR  results  are  used  to  investigate  the  diffusion  process  of 
a  chain  in  a  melt. 


i  '  — 


The  iotid-tike  &fU.n-&y&tejn  fiteponbe. 


i)  When  the  temperature  of  a  polymer  system  is  raised  starting  from 
its  glassy  state,  the  resonance  line-width  of  nuclear  spins  is  found  to  rapi¬ 
dly  decrease  through  the  glass  transition  temperature  ;  then  two  cases  must 
be  considered.  For  a  polymer  system  made  of  short  chains,  the  spin-system  made 
of  short  chains,  the  spin-system  has  a  liquid-like  response  (narrow  line  ana¬ 
logous  to  that  observed  on  conventional  liquids).  Whereas  for  a  polymer  system 
made  of  long  chains  the  spin-system  exhibits  a  solid-like  response  »  the  reso¬ 
nance  line  h,  5x1 02  Hz)  is  of  course  much  narrower  than  that  observed  in  the 
glassy  state.  This  is  the  first  basic  NMR  result.  The  solid-like  response  is 

easily  controlled  from  a  spectrum  narrowing  effect  induced  by  a  sample  rotation 

g 

around  an  axis  perpendicular  to  the  direction  of  the  steady  magnetic  field  . 
The  dL6e.nta.nqleme.irvt  NMR  tAamition  cuAve 

ii)  At  a  given  temperature,  a  transition  of  NMR  properties  from  a  solid¬ 
like  spin-system  response  to  a  liquid-like  one  is  induced  by  decreasing  tl  e 
chain  molecular  weight  j  correspondingly,  the  resonance  line-width  is  found  to 
decrease  from  about  5x102  Hz  to  about  10  Hz  according  to  a  reasonably  sharp 
curve  which  we  called  a  disentanglement  NMR  transition  curve  because  the  strong 

chain  length  dependence  of  the  resonance  line-width  is  necessarily  associated 

10 

inth  an  increase  of  the  rate  of  dissociation  of  entanglements  .  This  is  the 
second  basic  NMR  result^Fig 

III  -  RESIDUAL  DIPOLE-DIPOLE  INTERACTIONS 


Before  explaining  how  the  chain  migration  process  may  be  involved  in 
these  basic  NMR  properties,  it  is  worth  emphasizing  that  magnetic  relaxation 
processes  observed  on  most  polymer  systems  are  induced  by  dipole-dipole  inter¬ 
actions  of  nuclei  j  they  are  tensorial  functions  of  nucleus  coordinates  j  they 
strongly  depend  upon  the  distance  [r|  between  two  nuclei  and  also  upon  the 
orientation  of  the  vector  Joining  two  nuclei  with  respect  to  the  steady  magne¬ 
tic  field  direction.  For  the  sake  of  simplicity,  the  energy  of  dipole-dipole 
interaction  is  written  as  : 

W  «  (3  cos20  -  1 )/r3  (4) 

(8  is  one  of  the  angular' coordinates  of  r). 


* 
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Most  polymer  chains  are  made  of  proton  pairs  and  methyl  groups.  Within  a  pro¬ 
ton  pair  or  a  {CH3)  group,  the  distance  between  nuclei  is  a  constant  and  only 

* 

angular  properties  are  observed,  translational  diffusion  effects  between  dif¬ 
ferent  solvents  are  eliminated. 

We  consider  that  the  solid-like  response  reflects  a  residual  energy 
of  dipole-dipole  interactions  because  a  partial  spectrum  narrowing  effect,  only, 
is  induced  by  macromolecular  motions.  Although  the  chain  diffusional  motion 
is  isotropic,  it  is  like  observing  an  apparent  non-isotropic  motion  because 
the  diffusion  of  a  chain  over  its  own  dimension  is  achieved  within  a  time 
interval  longer  than  the  NMR  time  scale  of  measurements. 

IV  -  A  BASIC  ASSUMPTION  :  A  TWO-STEP  MOTIONAL  AVERAGING 
Two  dtipeAA-ionA 


The  only  basic  assumption  necessary  to  account  for  the  solid- like 

response  and  the  disentanglement  NMR  transition  curve  is  to  consider  that 

there  is  a  two-step  motional  averaging  of  spin-spin  interactions  »  this  is 

supposed  to  result  from  a  wide  chain  relaxation  spectrum  consisting  of  two 

11 

well  separated  parts  :  two  dispersions  8^  and  £22*  This  hypothesis  is  in 
agreement  with  general  viscoelastic  properties  and  with  the  dynamic  screening 
effect  :  to  the  cut  in  space  correlations  also  corresponds  a  cut  in  time 
correlations. 

is  called  the  transition  spectrum  ;  it  is  associated  with  local 
motions  of  short  segments  >  it  may  be  studied  form  the  nuclear  magnetic  spin- 
lattice  relaxation.  While  is  called  the  terminal  spectrum  ;  it  is  asso¬ 

ciated  with  the  collective  diffusion  of  all  parts  of  a  chain  ;  it  strongly 
depends  upon  the  chain  molecular  weight. 

Lou  o£  memory  o£  oAsLuttcvtioru 

The  memory  of  orientations  of  a  given  monomeric  unit  observed  through 
dipole-dipole  interactions, 3jfg,  is  lost  in  two  steps  i  the  memory  function  : 

♦  (t)  «  '<^(t)#0{o)> 


« 


(5) 


B.- 


i 


Is  splltted  In  two  parts,  with  : 


♦(o)  -  <(*£<»)  -  <3£>n  )2>„  *  |4): 


i)  There  is  a  fast  decay  of  the  memory  function  associated  with  the 
spectrum  i  this  leads  to  a  residual  energy  <  calculated  within  a  dynamic 

screening  domain,  l.e.  over  short  segment  motions  corresponding  to  the  0^  spec* 


trum(  Tig.  l)  . 


11)  There  is  a  long  decay  of  $Ct)  associated  with  the  spectrum. 


SpecX/ium  mAAouiina  conditions 


Two  conditions  must  be  fulfilled  to  observe  a  complete  spectrum  narro¬ 
wing  effect  : 

-  /tc  ^  1.  for  the  spectrum  [7] 


ii)  Tre  <  1.  for  the  spectrum 


xc  is  the  longest  relaxation  time  of  j  and  I  plays  the  role  of  an 

internal  NMR  reference  frequency  which  can  be  used  to  monitor  the  terminal 
relaxation  spectrum 


A  tempon.aA.tJ  network  stAuctuAZ 


Finally,  all  these  properties  may  be  pictured  by  considering  there 
exists  a  temporary  network  structure  made  of  temporary  submolecules.  A  sub- 
molecule  is  a  temporary  object  characterized  by  : 


1}  its  size  defined  from  the  dynamic  screening  length  O0  «  Ne’  , 


Ne  «  Ct1  rt° 


UTMWite 


7.- 


11)  Its  life-time  identified  with  the  relaxation  time  of  disentanglement 

TRbp  *  C"3  '  1,0> 

For  polymer  systems  made  of  long  chains,  we  only  observe  NMR  properties  of  the 
temporary  network  structure,  while  for  polymer  systems  made  of  short  chains 
we  observe  the  dynamical  dissociation  of  the  network  structure. 

Characteristic  NMR  properties  of  this  network  structure  must  be  now 
determined. 

V  -  RESIDUAL  V1P0LE-V1P0LE  ENERGY  UTTH1H  A  SUBHOLECULE 
Tampofiaiu  aveAage.  ofUtrUcutionaJt  omdv i 


Whatever  the  complex  nature  of  all  topological  constraints  exerted  on 
a  submolecule,  we  suppose  that  it  is  fully  characterized  at  any  time,  t,  by 
its  temporary  end-to-end  vector  rB(t)  j  it  is  supposed  to  be  the  only  rele¬ 
vant  thermodynamic  variable  describing  a  submolecule.  The  residual  energy  of 

dipole-dipole  interaction  corresponding  to  a  given  end-to-end  vector,  r„(t), 

12  13 

can  be  calculated  for  a  proton  pair  or  a  methyl  group  *  : 

E(rg)  -  I  «  <3  CO20  -  1>_  Cl  1 1 

1  re 

ECr  )  «  r  (  _1  2)r  /Ne*  M2) 

“  “  0 

The  non-zero  end-to-end  vector  rB(t),  induces  a  reduction  of  the  chain 

entropy  j  it  is  like  considering  that  every  monomeric  unit  experiences  an 

additional  potential  energy  hinderning  its  isomerisation  process.  This  entropy 

14  15 

reduction  effect  can  be  observed  on  polymeric  gels  ’  and  on  polymer  sys¬ 
tems  with  surface  interactions,  too. 

Formula  C12)  actually  describes  a  transfer  of  localisation  of  NMR  pro¬ 
perties  from  a  local  space  scale  to  a  semi-local  one  defined  by  re(t)  :  there 
is  a  temporary  average  orientational  order  of  monomeric  units.  This  order  is 
a  small  effect  since  it  reflects  a  strength  of  correlation  equal  to  about  102  Hz 
while  the  strongest  correlation  or  orientation  (in  a  glass)  is  about  10s  Hz. 

We  are  not  interested  anymore  in  the  memory  of  orientation  of  a  monomeric  unit 
but  in  the  loss  of  memory  of  orientation  of  rB(t)  vectors. 


1  /  H 


CaJUhKated  qvU  :  iupM.po6W.on  MR  ptopvUy 


NMR  properties  may  be  shown  to  obey  the  key  formula  (12)  by  observing 
the  progressive  swelling  process  of  elementary  chains  in  calibrated  covalent 
gels.  End-to-end  vectors  are  identified  with  vectors  joining  two  consecutive 
cross-link  points  j  they  hardly  vary  with  time  in  a  covalent  gel  corresponding 
to  a  permanent  network  structure.  The  characteristic  dependence  of  E(r0)  with 
respect  to  the  submolecule  (elementary  chain)  end-to-end  vector  re  is  observed 
by  swelling  the  covalent  gel.  Applying  a  packing  condition  to  partly  swollen 
elementary  chains.  |r0|  may  be  expressed  as  a  function  of  the  swelling  ratio, 
q.  according  to  the  formula  : 


(13) 


q  is  defined  as  the  ratio  of  the  volume  V  of  the  swollen  gel  over  the  volume  V 


of  the  dry  gel  (q  «  V/V0)^^«,  a). 


To  observe  the  characteristic  dependence  of  E  with  respect  to  the  num¬ 
ber  of  bonds  of  a  submolecule  (elementary  chain),  Ns  is  varied  by  changing  the 

15 

concentration  of  synthesis,  v  .  of  covalent  gels  : 

c 


N 


-5/4 

v_ 


(14) 


» 

Combining  (13)  and  (14).  the  residual  energy  of  dipole-dipole  interaction 
should  vary  as 


E  « 


(15) 


2/3  5/3 

Using  the  reduced  variable  q  v  ,  transverse  relaxation  rates  measured 

c 

on  calibrated  polydlmethylslloxane  gels  have  been  shown  to  obey  a  superposition 

15  ■» 

property  .  Several  concentrations  of  synthesis  v  ■  .84,  .74  and  .46  g/cmJ 

c 

and  two  swelling  agents  (toluene  and  chloroform)  were  used.  The  methyl  group 
is  well  appropriate  to  the  study  of  submolecule  properties  because  the  fast 
rotation  around  its  C-axls  eliminate  most  of  the  dipole-dipole  interaction.  The 
remaining  part  only  depends  upon  the  orientation  of  the  ?-axis  i.e*.  of  skeletal 


bonds. 
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Entangled  chain*  ;  a  Aoggjgg hJALon  pnopyity 

NMR  properties  may  also  be  shown  to  obey  formula  (12.)  by  Increasing 
the  mesh  size  of  the  temporary  network  structure.  The  number  of  bonds.  Ne.  In 
a  submolecule  is  increased  by  slightly  diluting  concentrated  solutions  of  long 
polymer  chains  according  to  formula  ( 9)  j  while  the  rate  of  dissociation  of 
entanglements  is  not  high  enough  to  induce  the  second  step  of  motional  avera¬ 
ging  process.  Therefore,  by  increasing  the  mesh  size  we  expect  to  only  observe 
a  reduction  of  the  residual  energy  of  spin-spin  interactions,  reflecting  more 
freedom  given  to  isomerisation  processes  of  monomeric  units.  The  relaxation 
function  of  the  transverse  magnetization  is  expected  to  keep  its  mathematical 
structure^  4^. 

Relaxation  functions  of  protons  recorded  on  five  concentrated  solu¬ 
tions  C.45  <  C  <  .9  g/cm^J  of  long  polyisobutylene  chains  CM  106)  in  carbon 
disulfide,  have  been  recently  shown  to  obey  a  superposition  property,  by 

A  0 

applying  a  suitable  shift  factor  to  the  time  scale  .  All  relaxation  functions 
correspond  to  a  solid-like  response  j  the  superposition  property  reflects  less 
average  order  of  monomeric  units  j  the  chain  entropy  is  Increased  within  sub¬ 
molecules. 

VI  -  CHAIN  DYNAMICS 

InteAncU  nefienence  frieguency 

We  now  have  a  quantitative  way  to  characterize  NMR  properties  of 
the  temporary  network  structure,  from  the  residual  energy  of  dipole-dipole 
interactions  : 

“  V2  (16) 

“1  * 

Ty  will  serve  as  an  internal  reference  frequency. 

A  submolecule  is  then  considered  as  an  ephemeral  element  which  will 
be  used  as  a  probe  to  analyse  long  range  fluctuations  in  a  chain.  We  suppose 
that  the  two-step  motional  averaging  process  still  applies  when  the  network 
structure  dissociates  quickly  enough  to  Induce  a  complete  spectrum  narrowing 
effect.  The  second  step  of  this  effect  is  closely  related  to  the  loss  of  memory 

of  orientation  not  of  a  monomeric  unit  but  of  a  end-to-end  vector  of  a  sub- 

« 

molecule. 


When  the  chain  molecular  weight  is  decreased,  the  mesh  size  is  kept 
constant  but  the  rate  of  dissociation  of  the  network  structure  is  increased. 
The  residual  energy  of  dipole-dipole  interactions  is  now  considered  as  a  time 
function,  E(re(t)),  obeying  the  obvious  condition  : 

E2(ra(tn  -  Tv‘2  (17] 

at  any  time  t.  The  description  of  the  network  at  any  time.  t.  The  description 
of  dynamical  fluctuations  of  the  network  structure  is  a  complex  many  body  pro¬ 
blem  which  has  been  solved  until  now  within  a  mean  field  approximation.  The 
dlssoclaprocess  of  the  network  structure  is  supposed  to  be  closely  reflected 
by  the  diffusion  motion  of  a  single  chain  in  dynamical  interactions  with  all 
surrounding  chains.  Therefore,  it  is  considered  that  the  loss  of  memory  of 
orientation  of  a  submolecule  end-to-end  vector  may  be  associated  with  the 
terminal  spectrum  of  a  chain.  Two  main  models  have  been  proposed  until  now  to 
describe  dlffusional  motions  of  a  chain  in  a  melt. 

Rou&e  model 


According  to  the  Rouse  model,  the  memory  of  position  of  a  polymer 
molecule  is  lost  in  each  of  its  points,  at  any  time.  This  model  is  built  from 
linear  thermodynamic  fluctuations  governed  by  a  free  energy  : 


«  3  kT  Tl  ,/He 

J 


(18) 


associated  with  the  end-to-end  vector  r  ,  of  the  whole  chain  is  a  sum  of  all 

a,  J 

elementary  free  energies  F  ,.  Collective  motions  of  submolecule  end-to-vectors 
-►  J 

r  ,  (t)  are  described  from  normal  modes  determining  the  Q9(f1)  spectrum.  Any 

J  R  ^ 

relaxation  time  is  defined  from 
P 


T 


R 

P 


(19) 


p 

T.  is  the  terminal  relaxation  time  i  normal  modes  have  an  uniform  statistical 
R 

weight  i  t1  is  predicted  to  vary  as  the  square  of  the  chain  molecular  weight. 
It  is  currently  considered  that  this  model  well  applies  to  short  chains  in  a 
melt,  without  any  entanglements  fj'rt  S). 
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R zptation  model 

The  other  model  proposed  to  describe  the  dlffuslonal  motion  of  a 
chain  In  a  melt  Is  founded  on  the  tube  concept  and  the  reptation  motion.  The 
tube  concept  was  first  Introduced  by  Edwards  to  describe  statistical  fluctua¬ 
tions  of  chains,  at  equilibrium.  The  reptation  motion  of  a  chain  in  its  surroun¬ 
ding  tube  was  proposed  by  Oe  Gennes.  According  to  the  reptation  model,  the 
memory  of  orientation  of  submolecules  is  lost  at  tube  ends  only.  The  memory 
of  orientation  is  Kept  in  central  parts  of  the  tube  until  these  are  reached 
by  one  of  the  ends  of  the  chain  moving  backwards  and  forwards  in  a  random  way. 
along  the  tube.  The  (^(M)  terminal  spectrum  consists  of  a  series  of  modes  cha¬ 
racterized  by  relaxation  times  : 

Tp°  “  W2  P-  1.3.5.... 

3 

The  terminal  relaxation  time  Tpep  was  predicted  to  vary  as  M  instead 
of  M2.  The  statistical  weight  of  each  mode  is  proportional  to  p-2  j  this  gives 
a  negligible  weight  to  all  modes,  except  for  the  first  one.  The  reptation  model 
is  in  reasonable  agreement  with  viscoelastic  properties^  6  ). 

CompaxliQn  with.  NMR  KUutti 


The  principle  of  the  NMR  observation  of  the  diffusion  process  of  a 
single  chain  in  a  melt  is  to  shift  the  spectrum  through  the  internal  refe- 

rence  frequency  T^  .  Starting  from  very  long  chains  ^  Is  shifted  towards  short 
values  of  chain  relaxation  times  by  decreasing  the  chain  molecular  weight  17 . 
Then,  the  relaxation  rate  of  the  transverse  componant  of  the  nuclear  magnetiza¬ 
tion  is  decreased  according  to  the  disentanglement  NMR  transition  curve  discus¬ 
sed  in  section  I. 

The  relaxation  function  of  the  transverse  magnetization  has  been 
calculated  according  to  the  formula  : 

CLtt)  ■  exp{-aT  2(  E  T2  [exp(-t  T  1  ) 
v  p^g  p,q  p,q 


♦  t  t'1  -  1])} 

p.q 


(20) 


I 


with  T_1  ■  t"1  ♦  t"1  (21) 

P.q  P  P 

a  is  a  constant  depending  upon  the  chain  diffusion  model  chosen  to 
13  17 

describe  NMR  properties  *  .  The  best  agreement  with  experimental  NMR  results 

is  obtained  by  combining  the  chain  molecular  weight  dependence  of  the  terminal 
relaxation  time  T„  given  by  the  reptation  model,  with  the  Q2(M)  relaxation 

•'®P  -JQ  A  g  4ia 

spectrum  given  by  the  Rouse  model  '  ’  ' .  Values  of  the  terminal  relaxation 

-2  -2 

times  obtained  from  previously  reported  NMR  measurements  were  5  x  10  ,  2  x  10 

-2 

and  8  x  10  sec  polydimethyl siloxane  (POMS),  polystyrene  (PS)  and  polylsobu- 

•3 

tylene  (PIB)  chains,  respectively  i  the  concentrations  were  1.,  .53  and  .47  g/an~ 

respectively  ;  while  the  molecular  weights  were  2.2  x  105.  2.5  x  10J  and 

2.25  x  105,  respectively.  These  relaxation  times  are  about  ten  times  longer 

than  those  usually  estimated  from  viscoelastic  measurements  and  defined  by  the 

ratio  n  /G°  »  n  is  the  zero  shear  rate  viscosity  and  G°  is  the  plateau  modulus. 

O  N  O  N 

13 

Results  observed  on  protons  were  found  to  be  similar  to  those  observed  on  C 
nuclei  in  natural  abundance.  Such  a  property  clearly  shows  that  magnetic  inter¬ 
actions  between  nuclei  located  on  different  chain  segments  have  negligible 
effects.  These  interactions  are  averaged  to  zero  by  local  motions  of  short 
segments  in  dynamic  screening  domains.  Therefore,  a  linear  distribution  nuclear 
magnetization  is  defined  along  polymer  chains  in  a  melt.  Consequently,  it  is 
possible  to  observe  dynamical  properties  of  a  single  chain  although  it  is  in 
dynamical  Interactions  with  all  other  surrounding  chains. 

CONCLUSION 

Long  range  fluctuations  of  polymer  chains  in  a  melt  observed  from 
NMR  correspond  to  the  loss  of  memory  of  orientation  of  temporary  submolecules 
associated  with  dynamic  screening  domains.  The  cirterion  of  isotropy  is  defined 
from  the  residual  energy  of  tensorlal  interactions  of  nuclear  spins  within 
screening  domains.  More  details  about  the  reptation  model  would  be  probably 
necessary  to  actually  account  for  all  W1R  results.  Also,  a  more  tedious  charac¬ 
terization  of  the  spir-syBtem  response  would  probably  lead  to  a  better  accuracy 
in  the  determination  of  the  chain  molecular  weight  dependence  of  terminal  rela¬ 
xation  times. 


x 

/ 


Finally,  it  is  worth  emphasizing  that  the  NMR  approach  to  the 
observation  of  a  chain  diffusion  process  necessarily  leads  to,  the  analysis  of 
the  whole  terminal  relaxation  spectrum  whefeas  measurements  of  a  diffusion 
constant  do  not  give  any  detailed  information  about  this  spectrum. 
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FIGURE  CAPTIONS 


A  schematic  disentanglement  NMR  transition  curve 

The  two-part  correlation  function  of  dipole-dipole  interactions, 
associated  with  the  two-dispersion  chain  relaxation  spectrum 

A  schematic  progressive  swelling  of  elementary  chains  in  a 
covalent  gel 

The  mesh  size  of  the  temporary  network  structure  is  Increased  by 
slightly  diluting  the  polymer  system  ;  the  number  of  bonds  Ne  in  a 
submolecule  is  increased 

According  to  the  Rouse  model,  the  memory  of  position  of  a  chain  is 
lost  in  each  of  its  points  at  any  time 

According  to  the  reptation  model  the  memory  orientation  is  loj£  at 
tube  ends. 

The  (12(11]  terminal  spectrum  is  shifted  through  the  internal  reference 
frequency  T^  by  decreasing  the  chain  molecular  weight. 
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Abstract 


h  ■  r  '  •; 

*  The  basis  of  *H  and  13C  high-resolution  NMR  investigations  of  solid 

C  .  . 

polymers  is  outlined.  The  l3C  NMR  spectra  of  solid  syndiotactic  and 
isotactic  polypropene  are  discussed  and  their  interpretation  in  terms  of 
conformation  and  chain-packing  effects  are  reviewed.  The  effects  of 
decreasing  temperature  on  the  13C  high-resolution  spectrum  of  an  annealed 
sample  of  isotactic  polypropene  is  described  and  interpreted  in  terms  of 
the  crystal  structure.  The  question  of  the  proportion  of  the  sample 
giving  rise  to  33C  signals  is  addressed  and  some  results  reported.  The 
main  cause  for  observing  only  part  of  the  total  sample  is  shown  to  be  the 

t  .  l 

rotating  frame  spin-lattice  relaxation  behaviour.  The  ,*H  spin-lattice 
relaxation  and  spectral  characteristics  of  a  number  of  polyolefine  samples 
are  summarised  and  the  role  of  spin-diffusion  discussed.  '  The 

a 

heterogeneity  revealed  by  the  multicomponent  Tlp  behaviour  has 
been  modeled  numerically  using  a  computer  and  the  conclusions  are 
outlined.  The  recovery  behaviour  of  poly-1 -butene,  following  conversion 
to  a  meta-stable  form  by  heating,  is  observed  using  both  cross-polarization 
and  single-pulse  excitation  techniques.  In  addition,  an  example  of  the 
investigation  of  a  propene/ethylene  copolymer  by  both  13C  and  methods  is 
described  to  Illustrate  the  interdependence  of  the  techniques. 
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1.  Introduction 

NMR  studies  of  solid  polymers  have  a  long  history  and  a  considerable 
volume  of  literature.  In  the  main,  until  about  a  decade  ago,  these 

investigations  relied  on  the  use  of  XH  broadline  spectra  and  associated 
measurements  of  spectral  moments  and  spin-lattice  relaxation  behaviour. 
More  recently,  other  nuclei  have  become  important,  largely  because  of 
advances  in  technology  and  developments  in  NMR  techniques.  In  particular 
the  use  of  2H  and  13C  has  become  widespread.  In  this  paper  we  describe 
some  work  carried  out  in  our  laboratory  which  involves  investigations  of 
solid  polyolefines.  It  is  presented  as  illustrative  of  some  of  the 
possibilities  of  both  13C  high-resolution  NMR  spectroscopy  and  *H 
broad-line  investigations.  No  attempt  is  made  at  a  review  of  other  work 
in  the  literature  and  the  interested  reader  is  referred  to  articles 
elsewhere  for  such  reviews  [1-4], 


2.  Experimental 


2.1  Samples 


(a)  Polyethylene:  Two  samples  of  high-density  PE  will  be  mentioned 
PEI  is  a  single-crystal  mat  obtained  by  slow  crystallisation  of  an 
0.1%  w/w  solution  in  xylene  at  70°C  followed  by  washing  with  methanol 
and  drying  under  vacuum  at  ambient  temperature.  The  second  sample, 
PEI I ,  is  a  melt-crystallised  sample,  annealed  at  403K  for  30  minutes. 

(b)  Polypropene:  The  Isotactic  polypropene  used  was  a  commercially 
available  grade,  manufactured  by  ICI  pic,  having  a  melt  flow  index  of 
20  and  an  isotactic  content  in  excess  of  97%  as  determined  by  proton 
NMR  in  solution  via  the  racemic  diad  concentration  (5).  PPI  was 
fabricated  as  a  thin  film  and  was  annealed  at  433K  for  60  minutes. 


PPII  was  also  a  thin  film  obtained  by  rapid  quenching 

from  the  melt  to  273K  followed  by  brief  heating  to  373K  to  remove 

any  of  the  smectic  form. 

(c)  Polybut-l-ene:  This  was  a  commercially  produced  sample  with 

an  X-ray  determined  crytallinity  0^30%.  PBI  was  used  as  obtained 
whilst  PBI I  was  made  by  heating  a  sample  of  PBI  to  just  above  the 
melting  temperature  (~380K)  followed  by  rapid  cooling  to  room 
temperature. 

(d)  A  propene/ethylene  copolymer:  a  block  copolymer  comprising 
85%  propene  as  a  main  block  and  an  end  block  made  with  a  50-50 
mixture  of  ethylene  and  propene.  The  polymer  was  examined  as  made 
(PECI)  but  was  also  subjected  to  extraction  with  boiling  heptane 
giving  two  further  samples,  the  heptane  soluble  (PECII)  and  insoluble 
(PECI I I )  fractions. 

2.2  NMR  measurements 

XH  NMR  measurements  were  made  at  two  operating  frequencies  of  60MHz 
and  200MHz  using  two  separate  spectrometers.  The  measurements  at  60MHz 
were  made  with  a  conventional  high-power  pulse  spectrometer  operating  with 
a  low-resolution  magnet.  Provision  was  available  for  adjustment  of  the 
r.f.  field  amplitude  which,  for  spin-locking  measurements,  was  generally  40 
or  60kHz  (=  -ylH  B^*).  Data  acquisition  was  by  means  of  a  digital 
signal  averager.  Samples  were  contained  in  10mm  o.d.  flat-bottomed  glass 
tubes  with  the  temperature  controlled  by  a  conventional  gas-flow  system, 
settable  to  +5°C. 

The  200MHz  *H  measurements  were  made  on  a  computer-controlled  double¬ 
resonance  spectrometer,  utilising  the  high-frequency  channel.  The 
spectrometer  is  based  around  an  Oxford  Instruments  wide-bore  4.7T  high- 
resolution  magnet  and  is  operated  by  a  Nicolet  1180  computer  system  and 
peripherals.  The  two  frequency  channels  each  allow  for  four  pulses  of 
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separately  adjustable  r.f.  phase  at  up  to  1KW  output  power,  the  latter 
being  under  computer  control.  Each  channel  has  a  broadband  quadrature 
receiver  followed  by  a  dual -channel ,  4-pole  Bessel -function  filter  of 
variable  bandwidth.  The  *H  spectra  were  recorded  using  a  digital  fast 
transient  recorder  operating  at  up  to  20MHz  with  8-bit  resolution. 

Normally  a  sweep  width  of  312.5  kHz  was  used  (3.2ys  dwell).  Samples  were 
contained  in  5mn  tubes  placed  in  the  horizontal  solenoid  coil  of  a  Bruker  Z 
32HP  probe.  Generally,  90°  pulse  lengths  of  the  order  of  lys  were 
employed  with  lower  r.f.  fields  being  used  for  long  spin-locking  pulses 
(40-60kHz). 

13C  high-resolution  spectra  were  measured  using  dipolar-decoupling 
(60kHz  decoupling  fields)  and  magic-angle  sample  rotation  (MAR)  (6).  The 
13C  signal  was  either  obtained  via  cross-polarization  from  protons  (7,8)  or 
by  normal  spin-lattice  polarization.  Some  of  the  13C  spectra  were 
measured  on  the  double-resonance  system  described  above,  (i.e.  v013C  = 
50MHz)  using  Andrew-Beams  style  of  rotors  fabricated  from  Macor  and  Delrin 
(polyoxymethylene).  The  latter  formed  the  rotor  base  and,  if  sufficiently 
far  removed  from  the  coil,  gives  little  or  no  signal.  Sample  volumes  were 
approximately  0.3cm3.  The  probe  is  a  Bruker  CPMAR  probe. 

Other  13C  spectra  were  obtained  on  a  spectrometer  operating  at 
90MHz/22.63  MHz  based  on  an  electromagnet.  The  probe  system  is  a  new 
design  based  on  the  spinning  system  described  elsewhere  (9,10).  It 
employed  cylindrical  sample  containers  made  of  glass  or  Macor  of  8imi  o.d. 
Only  the  sample  plus  coil  is  in  the  magnet  and  the  probe  is  remotely  tuned 
and  matched  for  both  frequencies.  Once  the  magic-angle  is  set  it  remains 
constant  even  on  changing  samples, uni  ike  the  Andrew-Beams  system  used  in 
our  Bruker  probe  where  it  is  necessary  to  include  some  solid  KBr  in  the 
sample  and  observe  the  79Br  resonance  to  allow  adjustment  of  the  angle  on 
each  sample  (11). 
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The  remotely-tuned  probe  on  the  90/22.63  MHz  spectrometer  allows  for 
variation  of  the  sample  temperature  whilst  performing  CPMAR  experiments  and 
some  results  are  given  later  using  this  facility. 

3.  General  Considerations 

Most  solid  polymers  are  structurally  heterogeneous  and  the 
semi-crystalline  polyolefins  are  no  exception  to  this.  This  heterogeneity 
reveals  itself  in  NMR  studies  of  these  materials  in  a  number  of  ways. 

Spectra  may  comprise  a  number  of  superimposed  lines  of  differing  width  and 
shape  (12-15)  and  spin-lattice  relaxation  may  be  more  complex  than  a  simple 
single  exponential  process.  These  LH  NMR  properties  are  usually  dominated 
by  the  magnetic  dipolar  couplings  between  the  protons  which  produce 
spectral  lineshapes  and  relaxation  behaviour  characteristic  of  the  spatial 
distributions  of  the  protons  and  their  relative  movements  due  to  thermal 
motions.  Thus  at  a  given  temperature,  crystalline  and  disordered  regions 
of  a  polymer  will  tend  to  give  rise  to  broader  and  narrower  resonances 
respectively  whilst  spin  lattice  relaxation  (T.  and  T.  ),  which 

i  ip 

depends  on  motional  frequencies  of  the  order  of  the  proton  precession 
frequency  in  B0  (typically  of  order  107-108  Hz)  and  BL  (typically  lO'MO5 
Hz)  respectively  for  T^  and  T^p,  can  vary  over  a  considerable  range 
of  values.  Generally,  the  rotating  frame  relaxation  is  much  more 
efficient  than  the  T^  process  in  solids.  All  the  above  Is  complicated 
by  the  fact  that  the  proton  dlplolar  coupling  produces  a  spatial  transport 
of  nuclear  magnetization,  called  spin-diffusion,  which  couples  together  the 
relaxation  behaviour  of  different  regions.  Thus,  as  is  shown  below,  a 
single  uniform  T^  behaviour  may  be  observed  when  T^pe  and 
spectra  Indicate  heterogeneity. 

Apart  from  the  Intrinsic  value  of  these  effects  in  the  study  of  solid 
polymers,  which  we  shall  give  examples  of  below,  they  are  important  in 


the  use  of  13C  high-resolution  NMR  since,  often,  the  work  done  with  this 
nucleus  involves  cross-polarization,  (CP),  in  which  the  13C  signal  is 
generated  from  the  XH  spin  system.  Two  typical  double-resonance  pulse 
sequences  used  for  13C  high-resolution  NMR  in  solids  are  shown  in  Figure 
1.  The  cross-polarization  sequence  (figure  la)  has  the  following 
features.  A  90x°  pulse  is  applied  on  the  channel  and  the 
magnetization  is  spin-locked  with  a  field  Bly.  The  delay  tj  is 
usually  zero  but  is  indicated  here  to  suggest  the  possibility  of  using  T2 
as  a  means  of  selectively  retaining  only  part  of  the  proton  magnetization 
for  subsequent  use.  A  further  delay,  t2,  allows  discrimination  or 
selection  by  means  of  differing  3H  T^p  values  before  a  contact  pulse  on 
the  13C  channel  produces  a  13C  signal  by  polarization  transfer  from  the 
protons.  It  is  clear  that  13C  signal  can  only  be  dervied  from  the  proton 
magnetization  remaining  after  both  the  T2  and  T-j  effects  in  times  tx 
and  t2  have  acted.  A  further  delay,  T3,  allows  for  discrimination  based 
on  13C-1H  dipolar  coupling  strengths  (16).  Following  CP  the  13C  signal 
is  recorded  in  the  presence  of  the  strong  dipolar  decoupling  field. 
Following  this,  when  t3  =  0,  a  90_x  pulse  flips  any  remaining 
magnetization  back  along  B0  (17).  The  time  TR  is  a  period  during 
which  the  system  undergoes  spin-lattice  relaxation  following  which  the 
procedure  is  repeated,  the  13C  signals  being  accumulated  in  the 
spectrometer  computer.  Figure  lb  illustrates  an  alternative  experiment  in 
which  the  13C  signal  Is  generated  via  its  own  spin-lattice  processes.  In 
this  case  the  time  TR  is  crucial  in  determining  which  13C  spins  give 
signals.  Only  those  for  which  T^^C)  TR  will  be  detected. 

These  sequences  are  usually  combined  with  magic  angle  sample  rotation  to 
yield  "liquid-like"  spectra  (6). 

With  these  general  comments  In  mind  we  now  discuss  a  number  of 
investigations  of  polyolefines  using  these  techniques. 
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4.  13C  spectra  of  polypropenes 

4.1  Effects  of  conformation 

The  13C  spectrum  of  solid  syndiotactic  polypropene  has  been  described 
elsewhere  (18).  The  interesting  feature  of  this  spectrum  is  a  1:1 
splitting  of  the  methylene  carbon  resonance  of  8.7ppm.  This  was 
interpreted  in  terms  of  the  conformation  of  this  molecule  in  the  solid, 
which  X-ray  diffraction  had  suggested  to  be  a  4*2/1  helix  in  which  there 
are  two  distinct  and  equally  probable  sites  for  methylene  groups:  one  on 
the  helix  axis,  the  other  on  either  side  of  this  axis.  Consideration  of 
shielding  effects  of  y-substituents  suggest  a  shift  difference  of~8ppm 
between  the  two  sites  which  is  close  to  the  observed  value.  In  addition 
the  higher  frequency  peak,  assigned  to  the  site  on  the  outside  of  the 
chain, shows  additional  broadening/splitting  consistent  with  interchain 
effects  (see  below). 

4.2  Effects  of  crystal  packing 

The  13C  spectra  of  samples  PPI  and  PPII  have  been  presented  elsewhere 
(19).  The  essential  features  of  the  spectra  of  these  isotactic 
polypropene  samples  are  that  PPI,  the  highly  annealed  material,  shows 
splittings  of  the  methyl  and  methylene  resonances  in  approximately  a  2:1 
ratio,  the  more  intense  peak  being  to  high-frequency  in  each  case.  Close 
examination  of  the  methine  resonance  shows  a  shoulder  in  the  same  sense. 

The  quenched  sample  does  not  show  these  splittings,  of  order  of  0.5-1.0ppm, 
and  the  spectra  Indicate  a  shift  of  resonance  Intensity  In  the  PPII 
methylene  and  methyl  regions  to  the  frequencies  characteristic  of  the  lower 
intensity  signals  in  these  regions  in  PPI.  The  interpretation  given  to 
these  observations  was  that  the  unit  cell  of  a  isotactic  polypropene 
contains  paired  left  and  right-handed  3:1  helices  and  that  this  generates 
distinguishable  sites  for  methyl,  methine  and  methylene  carbons  in  a 


2:1  ratio  by  virtue  of  inter-chain  interactions.  The  effect  of  quenching, 
on  this  basis,  would  seem  to  be  an  increase  in  resonance  intensity  at 
frequencies  characteristic  of  the  more  open  environment  associated  with 
interactions  between  chains  in  different  pairs  and  hence  the  suggestion 
that  quenching  traps  a  significant  fraction  of  the  chains  in  a  non-paired 
arrangement.  X-ray  diffraction  still  indicates  the  same  basic  structure  so 
these  effects  are  on  a  very  local  scale. 

4.3  Effects  of  decreasing  temperature 

Lyerla  et  al  have  reported  that  on  decreasing  the  temperature  of  solid 
isotactic  polypropene  the  13C  NMR  signal  from  the  methyl  group  broadens  due 
to  interference  of  the  modulation  of  the  13C-1H  dipolar  coupling  due  to  the 
thermal  motion  of  the  methyl  group  with  that  due  to  the  strong  proton 
decoupling  (20).  Figure  2  shows  the  temperature  variation  of  the  13C 
spectrum  of  the  annealed  sample,  PPI,  obtained  using  the  VT  probe  system 
mentined  in  the  experimental  section.  The  spectra  show  similar  behaviour 
to  that  reported  by  Lyerla  et  al  but,  owing  to  the  resolution  of  different 
peaks  in  this  sample,  various  other  features  are  apparent.  Firstly,  the 
methylene  signal  shows  little  change  over  the  range  of  temperatures 
covered, showing  the  maintenance  of  resolution  with  temperature  variations. 

Indeed,  the  Hartmann-Hahn  matching  and  the  general  electronic  properties 
of  the  probe  also  remained  unchanged.  Secondly,  it  is  clear  that,  as  the 
temperature  is  lowered  from  300K,  the  larger  of  the  two  methyl  resonances 
is  the  first  to  broaden.  This  can  be  taken  as  supporting  the  assignment 
of  this  splitting,  given  elsewhere  (19)  and  mentioned  above,  to  the 
existence  of  paired  helices  of  opposite  handedness.  The  methyls  within 
the  close-meshed  pair  of  helices  would  be  likely  to  experience  more 
hinderance  to  their  motion  than  those  on  the  outside,  this  leading  to  the 


observed  behaviour  with  temperature. 
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4.4  Quantitative  considerations 

Observations  of  the  13C  spectra  of  polybut-l-ene  led  us  to  determine 
the  proportion  of  the  total  carbon  in  the  polypropene  samples  that  was 
being  detected  by  the  usual  cross-polarization  experiment.  The  probe 
system  90°  pulse  length  was  calibrated  using  adamantane  and  the  efficiency 
of  the  cross-polarization  process  checked  using  hexamethyl benzene.  The 

proportion  of  carbon  was  then  determined  by  comparing  known  masses  of  PPI 
and  II  with  calibration  samples  of  hexamethyl benzene.  The  results 
obtained  indicated  that  for  PPI  and  PPII  approximately  80%  and  70% 
respectively  was  being  observed  (21).  The  reason  for  this  is  apparent 
when  the  XH  relaxation  properties  are  considered.  For  the  CP  process  the 
13C  signal  produced  depends  on  the  both  the  *H  Tjp  behaviour  and  the 
13C-1H  cross-relaxation  time  (22).  If  the  former  is  short  then  it  may  not 
be  possible  to  obtain  thermal  equilibrium  between  the  two  spin  systems 
before  the  proton  bath  is  significantly  heated  by  the  snin-lattice  process 
(Tn  ).  In  the  case  of  many  polymer  samples  lH  T.  behaviour  is 

Ip  IP 

multiexponential  in  character  and,  for  PPI  and  PPII  has  a  short-time 
component  which  Is  of  the  order  of  600-800ms  which,  in  turn,  is  associated 
with  a  narrow  (long  T2)  component  in  the  spectrum.  These  two  factors 
make  It  extremely  unlikely  that  any  significant  cross-polarization  arises 
from  that  part  of  the  proton  magnetization.  This  largely  accounts  for  the 
deficit  in  the  CP  experiments  mentioned  above. 

5.  *H  relaxation  and  spectra 

In  this  section  we  present  a  summary  of  results  which  have  been 
described  in  detail  elsewhere  (23).  The  proton  spin-lattice  relaxation  In 
each  of  PPI,  PEI I,  PPI,  PPII  and  PBI  at  room  temperature  is  described  by  a 
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single  characteristic  relaxation  time.  The  on-resonance  rota  ting-frame 
spin-lattice  relaxation  (T^90),  on  the  other  hand,  requires  a  minimum 
of  three  exponential  processes  for  Its  description  in  each  case.  Of 
particular  note  is  the  fact  that  the  long  time  behaviour  of  the  T,  90 

ip 

process  is  strongly  dependent  on  the  physical  history  of  the  samples. 

Thus  for  PEI  and  PEII  the  long-time  component  relaxation  times  are  60  and 
190ms  respectively,  whilst  for  PPI  and  PPII  they  are  109  and  30ms 
respecti vely. 

Further  measurements  of  the  spin-lattice  relaxation  behaviour  in  the 
off-resonance  rotating  frame,  with  the  magnetization  spin-locked  at  the 
"magic  angle"  of  54.7°  with  respect  to  the  B0  (Z)  axis,  showed  that  the 
relaxation  became  a  lot  slower,  whilst  remaining  three-component. 

A  number  of  possible  explanations  are  possible  for  this  overall 
relaxation  behaviour  but  the  most  likely  is  that  the  three  relaxation 
components  relate  to  three  physically  distinct  regions  in  the  polymers  in 
which  the  protons  have  different  NMR  properties.  This  heterogeneity  of 
structure  In  semi -crystal! ine  polymers  is  well-known  (2,3)  and  the  number 
of  spectral  components  used  to  fit  experimental  data  has  varied.  For 
example,  Bergmann  (13)  and  Smith  et  al  (15)  have  used  three  components  to 
fit  spectral  llneshapes  but  Bergmann  has  more  recently  chosen  to  describe 
the  data  in  terms  of  a  crystalline  component  plus  an  amorphous  component 
descrlbable  by  a  lineshape  involving  a  distribution  of  correlation  times 
(24).  Analysis  of  spectra  in  this  way  to  yield  information  on 
crystallinities  etc.  Is  quite  successful  but  subject  to  a  degree  of 
arbitrariness  in  choice  of  lineshape  functions. 

Spin-lattice  relaxation,  on  the  other  hand,  is  subject  to  the 
complication  of  spin-diffusion  (23).  This,  as  stated  above,  is  a  spatial 
transport  of  magnetization  via  homonuclear  dipolar  couplings  and  leads  to  a 
coupling  of  the  spin-lattice  relaxation  of  spatially  separated  regions. 
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This  can  be  invoked  to  explain  why  single  T-j's  are  observed  and 
multiple  Tip90's"  *n  t*1®  first  place,  T.j,  being  dependent  on 
fluctuations  at  o)0  and  2wo  (aio=yB0)  Is  often  considerably  longer  than 
Tip90’  which  depends  on  low-frequency  fluctuations,  of  order  (=yB1). 

In  addition  the  spin-diffusion  coefficient,  Ds>  is  proportional  to  the 
secular  dipolar  Interaction  strength  which,  under  on-resonance  spin-locking 
is  scaled  by  0.5.  Thus  T^  gives  a  longer  timescale  for  a  faster 
diffusion  process  to  take  effect. 

c  .  70 

The  observations  that  Tj  was  on  a  considerably 

longer  timescale  than  T-|p90  indicated  that  even  for  Tlp90. 
spin-diffusion  is  causing  a  partial  mixing  of  the  intrinsic  relaxation 
properties.  To  illustrate  this,  as  reported  elsewhere  (23),  we  have 
simulated  the  relaxation  behaviour  for  a  three  region  system  allowing  for 
spin-diffusion  in  one  dimension.  The  equation  governing  this  behaviour  is 

where  a  labels  the  magnetization  component  (e.g.  z  for  T^ ,  etc.),  q  the 

region  (q  =  1,2  and  3)  and  j  the  relaxation  process  (l,lp.  etc.). 

MzCq  =  Mo  otherwise  Meq  =  0.  The  spin-diffusion  coefficients 

and  relaxation  rates  are  different  In  each  region  (hence  the  label  q)  but 

constant  within  each  region  (not  a  function  of  x).  Numerical  solutions  of 

this  equation  using  values  for  the  sizes  of  crystalline  lamellae  etc.  taken 

from  other  sources  demonstrated  that  the  observed  behaviour  of  T.  90 

Ip 

could  be  explained  on  the  basis  of  the  model.  In  particular.  It 
demonstrated  that  the  long-time  Tj^90  components  largely  reflect  the 
time  for  magnetization  to  diffuse  out  of  the  crystalline  regions  Into  the 
less  ordered  regions  which,  for  these  materials  at  the  temperature  used, 
act  as  relaxation  sinks. 
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Other  experiments  involving  partially  relaxed  proton  spectra  have  been  i- 

t  carried  out  and,  in  particular,  experiments  using  T,  behaviour  to 

•  p 

spatially  label  the  sample  have  been  reported  (23).  These  Involve,  inter 
alia,  a  preparation  period  in  which  the  system  is  allowed  to  relax 
t  (Tlp90)  unt^  only  the  long-time  component  remains.  This 

magnetization  is  then  placed  along  B0  (z)  to  allow  maximum  spin-diffusion 
and  the  progress  of  this  mixing  is  examined  via  the  spectrum  obtained  with 
f  a  90°  pulse.  These  experiments  again  confirm  the  overall  behaviour  and 

heterogeneity  discussed  above. 
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6.  Isotactic  Polybut-l-ene 

The  13C  spectrum  of  PBI  obtained  using  CPMAR  is  shown  in  Figure  3. 

This  spectrum  shows  splittings  of  the  methylene  resonances,  probably 
arising  from  the  details  of  the  crystal  structure  as  with  the  isotactic 
polypropene.  The  XH  T^  behaviour  of  this  material  (23)  is  again 
three  component  with  the  longest  component  being  of  the  order  of  243ms  and 
the  short-time  component  (Tlp~500vs)  being  a  much  larger  proportion  of 
the  total  signal  than  for  the  PE  or  PP  samples.  After  heating  to  the  melt 
and  cooling,  the  sample  was  examined  by  both  the  CP  and  SPE  experiments  to 
record  13C  spectra  at  various  times  after  the  transformation  to  the  11:3 
helix  structure  which  this  treatment  brings  about.  Figure  4  Illustrates 
the  results.  It  should  also  be  noted  that  the  lH  Tj  immediately 
after  the  transformation  Is  close  to  being  a  single  exponential  with 
Tlp~  400ys.  This  fact  clearly  explains  why  the  CP  experiment  falls  to 
give  a  good  spectrum  immediately  following  the  transformation.  On  the 
other  hand,  the  SPE  experiment  also  falls  to  give  a  well -resolved  spectrum 
and,  since  the  *H,  T^  Is  250ms  the  13C  T^s  are  unlikely  to  be 
very  long  so  this  is  perhaps  unexpected.  The  likely  explanation  is  that 
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the  chains  are  executing  substantial  motions  which  give  rise  to  large 
linewidth  contributions  arising  possibly  from  the  decoupler/thermal  motion 
interfence  effect  described  earlier. 

As  the  recovery  of  the  material  towards  the  3:1  helix  form  proceeds 
the  CP  and  SPE  experiments  increasingly  give  Improved  quality  spectra  both 
in  terms  of  signal-to-noise  and  resolution.  It  is  Interesting  that  the 
best  spectra  are  obtained  after  some  two  weeks  of  annealing  but  that  the 
final  spectra,  taken  after  two  months,  seem  to  have  lost  resolution.  This 
is  possibly  because  the  material  which  recrystallises  first  gives  better 
crystals  and  that  the  eventual  apparent  decrease  in  resolution  arises  from 
the  slower  conversion  and  crystallisation  of  less  perfect  chains. 

7.  Propene-ethylene  co-polymers 

We  have  undertaken  studies  of  a  number  of  propene-ethylene  co-polymers 
and  the  results  of  one  such  investigation  are  given  here.  The  samples  were 
as  described  in  the  experimental  section.  Figure  5  shows  the  spectra 
(normal  and  partially  Tjp90  relaxed)  for  PECII  and  III,  the  heptane 
soluble  and  insoluble  fractions  respectively.  These  show  that  PECII I 
behaves  rather  similarly  to  PPI  and  II  as  described  elsewhere  (23). 

PECII,  however,  shows  the  presence  of  a  very  sharp  spectrum  with  a  small 
underlying  broad  compnent  which  Is  revealed  in  the  partially  relaxed 
spectrum.  Figure  6  Illustrates  the  the  13C  spectra  for  all  three  samples 
under  different  conditions.  Figure  6(a)  shows  the  standard  CPMAR 
experiment  for  PECI.  It  Is  essentially  the  spectrum  of  Isotactic 
polypropene  (19)  although  the  methylene  and  methyl  peaks  seem  more 
symmetric.  Use  of  the  delayed  contact  experiment  produces  changes  In  this 
spectrum  which  are  not  presented  here  but  which  indicate  a  composite  nature 
of  the  spectrum.  Figure  6(b)  shows  the  spectra  obtained  by  means 
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of  the  CPMAR  and  SPE  experiments  on  PEC III,  the  heptane  Insolubles.  These 
exhibit  some  distinctive  features.  The  methylene  resonance  in  the  CP 
spectrum  shows  some  signs  of  resolved  splitting  although  the  methyl  still 
does  not.  The  SPE  experiment,  with  a  recycle  time  of  lsec,  gives  a 
dominant,  unsplit  methyl  resonance  and  smaller  methine  and  methylenes. 

This  arises  because  the  methyl  groups  in  the  crystalline  portions  of  the 
main  polypropene  block  have  sufficient  motions  to  give  13C  T1  values 
short  enough  to  contribute  significantly  to  the  SPE  experiment  with  this 
recycle  time.  The  smaller  methine  and  methylene 'signals  are  probably 
associated  with  more  mobile  fractions  of  the  main  polypropene  block. 

Figure  6(c)  shows  the  CPMAR  and  SPE  spectra  for  PECII,  the  heptane 
solubles.  The  CPMAR  spectrum  is  essentially  that  of  isotactic  polypropene 
although  there  are  no  obvious  resolved  splittings  and  there  may  be 
additional  resonances  in  the  methylene  region.  The  SPE  spectrum  with  a 
recycle  time  of  lsec  however  shows  quite  a  different  pattern  of  lines. 

The  relative  intensities  of  the  lines  are  varied  by  changing  the  recycle 
period.  For  example  the  peak  at  the  frequency  corresponding  to  the 
methine  resonance  in  isotactic  polypropene  is  lost  when  the  recycle  time  is 
reduced  to  0.3sec.  The  detailed  interpretation  of  these  observations  will 
be  dealt  with  elsewhere. 
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Figure  2  The  variation  with  temperature  of  the  13C  high-resolution 

NMR  spectrum  of  annealed  isotactic  polypropene  (PPI).  Each 
spectrum  is  the  result  of  2000  FID's  acquired  with  a  60kHz 
decoupling  field,  5ms  contact  time  and  2s  recycle  time.  The 

spectrometer  frequencies  were  90MHz  for  *H  and  22.633MHz  for 
13C.  cijOuf ifctj,  Ikjt  J\AS  i-CC K 

Cyfcpb  Sp€tbnAj*%>  <P^oK 

Figure  3  The  13C  high-resolution  NMR  spectrum  of  solid  polybut-l-ene 

(PBI)  at  300K.  2000  FID's  were  added  and  were  obtained  with 

a  60kHz  decoupling  field  strength,  1ms  contact  and  Is  recycle 

time.  The  spectrometer  frequencies  were  as  for  Figure  2,‘fKe 

*e»k<KH:  $=(*•£>  ppm.,  C«3  j  S-Jt,- 8 frfcoW),  CH^a&dUi*): 

S=5I*,  CM ;  S=  SSi 

Figure  4  The  13C  NMR  spectra  of  solid  polybut-l-ene  (PBI I)  as  a  function 
of  time  following  conversion  from  PBI  by  heating/cooling  as 
described  in  the  text.  Single  pulse  excitation  (SPE) 
experiments  utilised  a  6.25ys  13C  90°  pulse  length  and  recycle 
time  of  2s.  All  experiments  employed  a  proton  decoupling  field 
strength  of  60kHz  and  the  CP  experiments  used  a  recycle  time 
time  of  1.2s.  The  spectra  are  Identified  as 
(a)  SPE,  Immediately  after  conversion;  (b)  CP,  Tc  *  0.6ms, 
t  =  4  hrs;  (c)  CP,  T  *  0.6ms,  t  *  8  hrs;  (d)  SPE,  t  *  12 

V 

hrs;  (e)  CP,  Tc  *  0.6ms,  t  *  16  hrs;  (f)  SPE,  t  ■  17  d*ys; 


Figure  5 


(g)  CP,  Tc  =  3.5  ms,  t  =  17  days;  (h)  CP,  Tc  =  3.5ms, 
t  =  50  days;  (i)  SPE,  t  =  50  days. 

NMR  spectra  of  (a)  PECII  and  (b)  PECIII,  the  heptane 
soluble  and  insoluble  fractions  of  the  propene/ethylene 
copolymer  PECI.  The  spectra  were  recorded  at  v0lH  = 

200MHz  following  spin-locking  *[(u>i/2ir)  =  40kHz^  for  a  time 
T  as  indicated.  Typically,  400  FID's  were  accumulated  with 
a  recycle  time  of  3s  and  a  solid  echo  was  used  to  measure  the 
signal  after  the  end  of  the  spin-lock. 

Figure  6  (a)  The  CPMAR  13C  high-resolution  NMR  spectrum  of  PECI. 

4000  FID's  were  accumulated  using  a  contact  time  of  1ms,  a 
recycle  time  of  Is  and  a  60kHz  decoupling  field  strength. 

(b)  The  13C  spectra  of  PECIII  obtained  using  SPE  (lower  trace) 
and  CP  (upper  trace).  The  recycle  time  for  both  experiments 
was  Is  with  other  conditions  as  for  (a). 

(c)  As  for  (b)  except  the  sample  is  PECII,  the  heptane  soluble 
fraction  of  PECI. 
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Summary 

- -  0 

Tt-rt.m.r\N  is  introduced  as  a  tool  for  investigating  slo w  molecular  motion 

in  the  glass  transition  region  of  amorphous  polymers.  In  particular,  we 
2  r,? 

comoare  spin  alignment  echo  spectra  of  chain  deuterated  polystyrene 
with  model  calculations  for  restricted  rotational  Brownian  motion.  Molecular 
motion  in  the  polystyrene-toluene  system  has  been  investigated  by  analysing 

p 

H  n.m.r.  of  partially  deuterated  polystyrene  and  toluene,  respectively. 

The  diluent  mobility  in  the  mixed  glass  has  been  decomposed  into  "solid" 
and  "liquid"  components  where  the  respective  average  correlation  times  differ 
by  more  than  5  decades.- 
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INTRODUCTION 


Previous  n.m.r.  investigations  of  the  glass  transition  in  polymers1  have 

orovided  information  on  relatively  raoid  motions  with  correlation  times 
-4 

below  10  s  at  temperatures  well  above  the  static  glass  transition 

2 

temperature  T  .  Recent  advances  of  H  n.m.r.  have  extended  the  dynamical 
range  to  slower  motions,  and  in  addition,  they  yield  specific  information 

Of 

uoon  the  type  of  the  molecular  reorientation  processes  involved.  In  two 
component  glasses,  there  is  the  further  opportunity  of  investigating  se¬ 
parately  the  molecular  dynamics  of  both  components  by  appliying  selective 
deuteration.  Atactic  polystyrene  has  for  decades  played  the  role  of  a 

"typical  example"  for  the  glass  transition  in  amorphous  polymers.  Thus,  we 
2 

are  also  studying  H  n.m.r.  in  polystyrene  samples  deuterated  at  the  chain 

(PS-d3)  and  at  the  phenyl  rings  (PS-d&),  respectively^  PS  and  toluene  (TOL) 

form  a  mixed  glass  over  the  whole  concentration  range  which  has  already  been 

investigated7  by  n.m.r.,  dielectric  relaxation,  and  thermal  methods 

(DSC,  DTA).  It  is  apparent  from  this  work  that  the  TOL  molecules  are  highly 

2 

mobile  at  temperatures  well  below  T^  of  the  mixture.  A  H  n.m.r.  study  where 
either  toluene  (TOL-d^,  TOL-d^)  or  polystyrene  (PS-d^,  PS-d^ )  are  deuterated 
should  provide  many  details  of  the  motion  in  the  glass  transition  region. 


In  the  following  sections  we  first  give  an  introduction  into  the  relation 
2 

of  H  n.m.r.  observables  with  molecular  reorientation,  for  details  we  refer 
2 

to  a  recent  review.  Furthermore,  we  attempt  to  show  the  extent  of  new  in- 

2 

formation  that  can  be  gained  from  H  n.m.r.  with  respect  to  the  glass 


transitidn,  taking  j|S  a>  a and  we  give  some  details  of  the 

1  gtr~~  ■ 

polystyrene -toluene  system.  ’ 
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2U 

H  n.m.r. 

In  a  rigid  solid  ,  the  coupling  of  the  deuteron  quadrupole  moment  with  the 

2 

axially  symmetric  field  gradient  of  a  C  -  H  bond  gives  rise  to  a  splitting 
2 

2  wq  of  the  H  Larmor  frequency  where 

2 

“Q  =  - gf 53  <3  cos2  8  "  (!) 

and  e  is  the  angle  of  the  C  -  H  bond  with  respect  to  the  external  magnetic 
field.  The  weighted  superposition  of  these  doublets  due  to  the  random 

orientations  in  the  glass  yields  the  characteristic  Pake  line  shape  of  the 

2  2  ? 

H  spectrum  having  a  width  of  ~  250  kHz  for  C  -  H  bonds.  Reorientation 

of  the  C  -  H  bond  can  be  described  by  a  stochastic  process 

fi(t)  =  cos  8 ( t)  (2) 

Q 

which  can  be  modelled  by  assuming  rate  equations 

P^/Rj.t)  =  Z  P(Qi/Qk,t)ity  (3) 

where  the  P(n-/R-,t)  are  conditional  probabilities  and  the  are  rates 

•  J  ■  J 

of  the  transitions  n.  ->•  n.  .  The  latter  can  be  chosen  in  accordance  with 

.  Q 

particular  assumptions  upon  the  reorientation  process.  It  can  be  shown 
that  the  line  shape  of  the  n.m.r.  spectrum  is  given  by  the  Fourier  transform 
of  the  free  induction  decay  (FID)  which  can  be  expressed  as 

S0(t)  =<(exp  i  /fc  u>Q(f)dt’)>  (4) 

where  the  average  can  be  evaluated  by  solving  a  rate  equation  similar  to 

equation  ,( 3}  containing  the  same  transition  rates  ir. .  of  the  model  for 

reorientation.  Equation  (4)  applies  in  the  whole  "slow  motion"  region  from 

the  rigid  Pake  SDectrum  to  the  raDid  motion  limit.  Detectable  line  shape 
2 

changes  of  H  spectra  occur  in  a  correlation  time  range  of  ■»»  1  -  20  us  . 
This  time  scale  can  be  extended  to  ~  200  ps  by  investigating  the  solid 
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echo  spectrum  which  can  be  formulated  as  the  Fourier  transform  of 

Sl(trTl)  =  <(expfi  -r  -  i  f  1wQ(t,)dt'  -  /tLQ(t')dt'}^  .  (5) 

is  the  distance  between  the  two  pulses  of  the  solid  echo  sequence  and 

tj  is  the  time  starting  from  the  echo  maximum  at  2t^  (see  Fig.  1). 

Equation  (5)  can  be  evaluated4  from  the  same  rate  equation  that  is  used  for 

solving  equation  (4).  As  a  matter  of  fact,  and  SQ  become  equal  in  the 

rigid  solid  limit  since  the  contributions  of  the  first  two  integrals  of 

equation  (5)  cancel  in  this  case.  However,  for  correlation  times  t  <  t, 

c  ~  1 

no  cancellation  occurs,  and  may  differ  considerably4  from  SQ.  The 
dynamical  range  of  the  solid  echo  technique  is  limited  by  the  transverse 
relaxation  time  T2  since  no  solid  echo  is  detectable  for  >>  T2.  In  most 
solid  polymers  T?  ~  200  -  500  us  for  2H  spectra.  It  should  be  noted  that 
T2  is  related  with  the  reciprocal  homogeneous  width  that  would  be  seen  in 
a  single  crystal,  c.f.  equation  (1).  Ultra  slow  motions  with  correlation  times 
up  to  w  I  s  can  be  investigated  by  application  of  the  Jeener-Broekaert  pulse 
sequence  (Fig.  1)  that  yields  a  spin  alignment  echo  after  the  3rd  Dulse.3 

O  C 

Alignment  echo  spectra  are  obtained  as  the  Fourier  transform  of  ’ 
s2(4,t1,t2)  =  (sin  [^(0)^]  sin  [Uq(t2)(t1  +  t£)]^ 

=  (cos  |(Ug(o)t2  +  £ug(o)  -  wq(t2)]  -  (  FID/*  .  (6) 

The  term^  ^FID^  is  of  minor  importance  in  the  context  of  the  present 
paper,  it  is  assumed  that  the  mixing  period  t2  is  much  longer  than  the 
evolution  and  detection  speriods,  respectively,  and.  that  n  changes  only  on 
the  time  scale  of  tj.  Thus,  t2  »  tj,t2  and  only  the  constant  orientations 
n(o)  and  £J(t2)  enter  in  equation  (6)  where  products  of  the  form  oigt  have 
replaced  the  integrals  ('ug(t,)dt'  of  equations  (4)  and  (5).  In  the  rigid 
solid  limit,  the  alignment  echo  S2(t2,Tj,r2)  is  a  trivial  superposition  of 


*■ 


! 


-  5  - 

a  solid  echo  and  the  term  <^F I D^>  .  Molecular  motion  is  seen  through  the 

difference 

Auq  =  0Jq ( o )  -  uq(t2)  =  (9  e2Qq/3fy  (cos2  eQ  -  cos2  9)  (7) 

It  is  remarkable  that  reorientation  by  small  sten  rotational  Brownian 

motion  causes  dramatic  changes  in  the  al  ignment echo  spectrum  provided  t2 

is  of  the  order  of  the  time  is  between  2  angular  steps  and  AwqTj  /v  it. 

Since  ts  is  much  smaller  than  the  rotational  correlation  time  tc  the 

alignment  echo  method  provides  an  unique  ODDortunity  for  distinguishing 

random  small  angular  step  reorientation  from  large  angular  jump  motion 

where  t  and  t  are  of  the  same  order  of  magnitude.  It  should  be  noted 
sc 

that  the  average  ^•••/>  of  equation  (6)  can  be  evaluated  from  numerical 
solutions  of  equation  (3)  where  small  step  rotational  Brownian  motion  is 
described  by  a  finite  approximation  of  the  diffusion  equation. 


BULK  POLYSTYRENE 

Our  extensive  study  of  partially  deuterated  PS-d^  and  PS-d^  has  already 
provided  a  wealth  of  information  upon  molecular  motions  above  and  below  T 

g 

2  10 

part  of  which  has  been  summarized  in  previous  reviews.  ’  In  this  section, 
we  whish  to  concentrate  on  the  nature  of  chain  motion  in  the  glass  transi¬ 
tion  region  where  new  aspects  have  evolved  from  an  analysis  of  alignment  echo 

spectra.  The  model  calculations  shown  in  Fig.  2  confirm  the  large  changes  at 

2 

times  t 2  <<  tc  where  the  C  -  H  bonds  can  only  reorient  within  a  small 
angular  range  of  a  few  degrees.  Furthermore,  they  reveal  a  surprising 
periodicity  for  t2  >  rc  with  a  peak  to  peak,  distance  proportional  to  tj1. 

Of  course,  the  assumption  of  constant  p(o)  during  the  evolution  period  x^ 
is  critical  for  the  appearance  of  this  spectrum.  On  the  other  hand,  the 
periodicity  can  be  suppressed  by  restricting  the  rotational  Brownian  motion 
to  a  limited  angular  range.  Thus,  the  alignment  echo  spectrum  for  rapid 
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restricted  diffusion  within  -  6°  (x«  >  5  t  )  is  similar  to  that  of  un- 
restricted  diffusion  for  x^  -  0.005  rc  where  the  angular  range  is  "re¬ 
stricted"  by  the  time  window  x^  of  the  experiment.  The  experimental  spectra 
of  Fig.  3  show  the  typical  behaviour  of  small  angle  reorientation  where 
at  380  K  the  angular  range  of  about  -  10°  is  attained  at  x^  =  40  ms. 

The  absence  of  periodicities  for  longer  values  rules  out  unrestricted 
rotational  diffusion.  Furthermore,  the  integral  of  the  theoretical  un¬ 
restricted  diffusion  spectra  essentially  vanishes  at  t^/t  >  0.1  since 
the  positive  and  negative  contributions  to  the  spectrum  almost  cancel.  The 
corresponding  experimental  intensities  decay  much  slower  for  long  values 
though  they  do  not  become  constant  as  predicted  from  the  model  of  restricted 
rotational  diffusion.  The  decay  of  the  alignment  echo  is  not  yet  fully  clari¬ 
fied,  however,  we  suggest  tentatively  that  only  a  fraction  of  the  C  -  H  bonds 
can  fully  reorient,  most  likely  through  conformational  changes, on  a  time 

scale  of  ~  1  s  and  this  fraction  decreases  rapidly  below  Tg.  At  higher 

2  2  10 

temperatures  T  >  400  K,  we  know  from  the  H  c.w.  absorption  spectra  ’ 
that  full  reorientation  of  the  C  -  H  bonds  occurs  by  essentially  a  single 
process  where  any  "distribution  of  correlation  times"  must  be  restricted 
to  less  than  one  decade.  Now,  it  turns  out  that  the  correlation  times  de¬ 
term  i red  from  the  slow  motion  line  shapes  are  at  least  an  order  of  magnitude 

larger  than  those  obtained  from  spin  lattice  relaxation  times  T^  in  this 
6  8 

temperature  range.  Apparently,  the  latter  are  determined  by  the  fast  re¬ 
stricted  rotational  diffusion  discussed  above.  Thus,  we  conclude  that  chain 
motion  in  the  glass  transition  region  can  be  described  by  a  rapid  rotational 

diffusion  process  wfjich  is  restricted  by  an  angular  range  of  approximately 

-  | 

-  10°  at  Tg  and  decreases  at  lower  temperatures.  Full  reorientation  of  the 

2 

C  -  H  bonds  occurs  by  a  slower  process.  It  should  be  noted  that  H  n.m.r. 
is  presently  the  only  experimental  method  for  detecting  the  reorientation 
of  a  single  bond  vector  in  the  glass  transition  region.  Photon  correlation 
spectra  as  well  as  dielectric  and  mechanical  relaxation  are  related  with 
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fluctuations  or  larger  volume  elements  where  the  cooperati vity  caused 
by  free  volume  redistributions  may  be  the  reason  for  finding  broad 
distributions  of  correlation  times. 

THE  POLYSTYRENE-TOLUENE  SYSTEM 
13 

It  is  well  known  that  the  glass  transition  T^  of  a  polymer  can  be 
considerably  reduced  by  adding  relatively  small  amounts  of  a  low  molecular 
weight  diluent.  If  phase  separation  occurs  due  to  thermodynamical  instabi¬ 
lity,  no  further  reduction  of  Tg  below  the  demixing  temperature  is  possible. 

However,  polystyrene  and  toluene  remain  miscible  over  the  whole  concentra- 

7  2 

tion  range  down  to  the  glass  transition  of  toluene.  Our  H  n.m.r.  study 
of  this  system  has  provided  information  upon  molecular  motion  that  may  be 
typical  for  any  glassy  mixture  of  a  polymer  and  a  small  molecule  diluent. 

In  bulk  PS-dg,  we  have  discovered  that  about  20%  of  the  phenyl  groups  per¬ 
form  180°  jumps  around  their  axis  to  the  backbone  even  at  temperatures 
2 

well  below  T  .  By  adding  toluene,  the  glass  transition  can  be  reduced  to 

low  temperatures.  We  find  that  the  fraction  of  flipping  phenyl  groups  becomes 

smaller  and  more  difficult  to  detect  for  diluent  concentrations  where  the 

solid  glass  regime  is  shifted  to  lower  temperatures.  There  is  no  indication 

of  large  amplitude  phenyl  group  libration  as  was  concluded  from  *H  wide  line 

7  14 

n.m.r.  by  Adachi  et  al .  ’  Thus,  we  could  use  PS-dg  for  investigating  the 
chain  motion  in  the  PS-TOL  system.  In  Fig.  4,  we  have  shown  the  T 2  values 
of  PS-dg  determined  by  measuring  the  echo  height  as  a  function  of  the  di¬ 
stance  Tj  between  the  two  pulses  of  the  solid  echo  sequence.  The  plot  of 
log  T2  versus  1/T  is  linear  at  temperatures  above  the  glass  transition,  and 
it  changes  to  an  almost  constant  value  at  a  temperature  T  p  close  to  the 
static  Tg  as  determined  from  other  experimental  methods/  cf.  Fig.  7.  In 
the  toluene  rich  region,  T2  decreases  as  the  temperature  is  lowered  below 
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T  D,  attains  a  minimum,  and  increases  again  to  the  rigid  solid  value. 

9  »■ 

Thus  a  second  temoerature  T'  D  can  be  estimated  as  is  shown  in  Fig.  4 

q  ,p  3 

for  a  solution  of  35.1%  PS-d^  in  T0L-dQ.  Apoarently,  there  is  a  temperature 
range  between  T  p  and  T'  p  where  the  chains  are  not  fully  immobilized 
due  to  the  abundance  of  highly  mobile  solvent  molecules  (see  below).  The 
solid  echo  line  shapes  of  PS-d^  will  be  discussed  in  a  later  publication.^ 

The  diluent  mobility  has  been  investigated  by  using  TOL-d^  or  TOL-d,..  No 

substantial  differences  between  these  toluenes  was  found  for  the  overall 

TOL  motion.  However,  the  rapid  methyl  group  rotation  in  TOL-d^  persists 

even  at  the  lowest  temperatures  of  our  experiments  thus  reducing  the  widths 
2 

of  all  H  spectra  to  one  third  of  the  corresponding  TOL-d^  values.  On  cooling 
a  TOL-dj  (or  TOL-dg)  solution  of  PS  below  the  temperature  T^  ^  shown  in 
Figs.  6  and  7,  a  solid  echo  can  be  detected  on  top  of  the  "liquid"  FID.  If 
this  signal  is  Fourier  transformed  with  respect  to  tj  starting  at  t  =  2  t.j 
(cf.  Fig.  1)  we  obtain  a  spectrum  as  shown  in  Fig.  5a.  The  central  portion 
originates  from  part  of  the  "liquid"  FID  that  starts  at  t  =  0.  It  should 
be  noted  that  the  second  90°  pulse  of  the  solid  echo  sequence  has  no  in¬ 
fluence  upon  the  FID  of  a  liquid  since  it  is  shifted  in  pase  by  90°  with 
respect  to  the  first  pulse.  In  Fig.  5b,  ij  has  been  increased  to  a  value 
where  the  FID  has  decayed  to  zero  (2  >>  T2  pID).  The  solid  echo  spectrum 

of  Fig.  5b  has  essentially  the  rigid  limit  Pake  line  shape.  Thus,  the  ro¬ 
tational  correlation  time  of  the  TOl-d^  molecules  that  give  rise  to  this 
spectrum  must  be  larger  than  100  us.  Since  we  can  detect  a  spin  alignment 
echo  after  50  ms,  the  fcorrelation  t fines  ire  gvet»  at>oye slow  time 
scale.  On  the  other  hand,  the  TOL-d^  molecules  contributing  to  the  "liquid" 
FID  have  correlation  times  below  ~1  ys.  Though  this  corresponds  to  re¬ 
orientation  in  a  rather  viscous  liquid  the  correlation  time  difference  of 
more  than  5  decades  provides  sufficient  evidence  for  qualitatively  different 
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behaviour.  We  have  tried  to  estimate  the  intensities  of  the  "liquid" 
and  "solid"  contributions  by  extrapolating  the  FID  to  t  +  0  and  the 
amplitude  of  the  solid  echo  to  -*■  0.  Since  the  solid  echo  amplitude 
become  non-exponential  at  short  times^  this  estimate  is  admittedly  crude. 

We  should  also  emphasize  that  the  decomposition  implies  no  statement  upon 
Dossible  correlation  time  distributions  within  the  "solid"  and  "liquid" 
components.  Nevertheless,  it  allows  for  some  important  conclusions  upon 
diluent  motion  in  the  mixed  glass.  The  TOL  solid  echo  becomes  observable 
above  T  D  at  temperatures  where  PS-dr  also  gives  rise  to  a  solid  echo 
(cf.  Fig.  4).  Thus,  the  mobility  of  the  "rigid"  TOL  molecules  is  linked 
to  that  of  the  PS  chains  at  T  <  T  cc.  There  are  even  indications  that 
it  is  slower  as  will  be  discussed  in  ref.  6.  On  lowering  the  temperature 
below  Tg  p,  increasing  amounts  of  TOL  "freeze"  towards  the  chains  until 
the  whole  mixture  becomes  a  solid  glass  at  Tg  pjp.  It  is  apparent  from 
Fig.  6  that  the  intensities  of  the  "solid"  and  "liquid"  components  add  to 
an  approximately  constant  total  intensity.  The  latter  has  been  scaled  in 
Fig.  6  by  an  arbitrary  factor  appropriate  for  presentation  of  the  data 
in  one  figure,  and  by  the  Boltzmann  factor  determining  the  trivial  tempe¬ 
rature  dependence  of  the  soin  populations.  Thus,  the  intensities  given  in 
Fig.  6  are  a  measure  of  the  "solid"  and  "liquid"  fractions  in  the  tempera¬ 
ture  region  Tg,FID  -  T  -  Tg,SE'  This  region  becomes  very  broad  for  large 
PS  concentrations.  Furthermore,  the  solid  echo  line  shapes  of  the  "solid" 

TOL  fraction  become  motional Iy_.iw rowed  for  PS  concentrations  above  ~  80% 

'*  ’  -  '  *  -  *  "='i  ?  •  -  :  - •  -  g 

indicating  the  abundance  of  "free  volume"  between  the  rigid  PS  chains. 

It  is  clearly  visible  in  Fig.  7  that  the  static  glass  transition  T^  (dotted 
line)  is  close  to  T  p  in  the  PS  rich  regime  indicating  that  the  solid  glass 
behaviour  is  governed  by  the  polymer.  On  the  other  hand,  the  TOL  glass  transi 
tion  takes  over  in  the  TOL  rich  regime  with  less  than  ~  50%  PS.  Here,  the 
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static  Tg  is  close  to  Tg  p  and  Tg  p and  some  chain  mobility  persists  in 

the  temoerature  region  T  D  >  T  >  T'  D.6  The  temperatures  T  n  and  T'  n 

g,P  -  g,r  r  g,P  g,P 

converge  at  very  high  TOL  concentrations  where  the  TOL  glass  transition 
could  be  determined  from  measurements  in  T0L-d3  having  a  similar  tem¬ 
perature  dependence  as  shown  in  Fig.  4  for  PS-d^.® 

Finally,  we  wish  to  note  that  different  diluent  mobilities  in  PS  have  also 

been  found  for  the  "dual  mode  sorption"  of  ammonia  as  investigated  by 
1  15 

H  n.m.r.  and  for  di-n-alkyl-phthalates  where  dielectric  relaxation  data 
were  interpreted  in  terms  of  three  different  a  processes.^ 


CONCLUSIONS 

2 

We  conclude  from  a  comparison  of  H  spin  alignment  echo  spectra  with  model 
calculations  that  chain  motion  in  the  glass  transition  region  of  PS  can  be 
described  by  a  relatively  rapid  rotational  diffusion  process  which  s  re¬ 
stricted  to  an  angular  range  of  approximately  -  10°  close  to  T  .  Full  re¬ 
orientation  by  a  slower  process  tends  to  become  motionally  heterogeneous  as 
the  glass  transition  is  approached  from  above. 

2 

H  n.m.r.  spectra  of  deuterated  toluene  in  the  PS-TOL  glass  can  be  inter¬ 
preted  in  terms  of  a  two  component  decomposition.  The  motion  of  the  "solid" 

component  is  linked  to  that  of  the  PS  chains  at  temperatures  below  a 
temperature  Tg  ^  where  a  solid  echo  becomes  observable.  On  lowering  the 
temperature,  the  fraction  of  "solid"  TOL  increases  until  the  whole  mixture 
becomes  a  solid  glass  a  temperature  Tg  pjp.  There  is^.a  rajtge  ...... 

Tg  FID  -  T  -  Tg  in  wliere  a  "liquid"  fraction  of  highly  mobile  TOL  mole¬ 

cules  exists  at  temperatures  below  the  "static"  glass  transition  temperature 
T  as  determined  by  dielectric  relaxation  or  thermal  methods. ^  T  is  close 

9  g 

to  T  ££  in  PS  rich  systems  and  close  to  T  in  systems  with  PS  concen- 


trations  below  50%  (see  Fig.  7).  In  the  TOL  rich  regime  we  can  determine 
a  temperature  Tg  p  >  Tg  where  the  chain  motion  becomes  highly  restricted 
and  a  "solid"  TOL  fraction  is  observable. 
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Figure  Captions 

Figure  1.  Notation  for  the  solid  echo,  and  the  generalized  Jeener-Broekaert 

3  pulse  sequence.  0  ,  <j>  ,  and  (|»  are  the  flip  angles  of  the  respective  h.f. 
y  ^ 

pulses. 

2 

Figure  2.  Theoretical  H  alignment  echo  spectra  for  the  model  of  rotational 
Brownian  motion. 

First  row:  unrestricted  diffusion  for  different  ratios  of  the  pulse 
distance  (Fig.  1)  over  the  rotational  correlation  time  t  . 

Second  row:  rapid  motion  limit  >  5)  for  diffusion  restricted  to 

the  angular  ranges  shown  in  the  figure. 

2 

Figure  3,  Experimental  H  alignment  echo  spectra  of  chain  deuterated 
polystyrene  above  T  . 

Figure  4.  Spin  relaxation  times  of  PS-dg  in  PS-TOL  mixed  glasses. 

Numbers  in  the  figure  denote  wt.%  PS-dg.  The  arrows  indicate  glass  transition 
temperatures  Tg  p  and  T^  p  (see  text). 

2 

Figure  5.  H  solid  echo  (a,b)  and  alignment  echo  (c)  spectra  of  TOL-d^  in  a 
mixed  glass  with  28.7  %  PS. 

2 

Figure  6.  H  intensities  (I  in  arbitrary  units,  see  text)  of  T0L-d3  (tri¬ 
angles)  and  TOL-dg  (circles,  squares)  in  mixed  glass  with  PS.  Numbers  tn  . 

the  figure  denote  wt.%  PS.  Open  and  full  symbols  denote  "liquid”  and 
"solid"  components,  respectively  (see  text). 
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Figure  7.  Glass  transition  temperatures  in  the  PS-TOL  system. 

Abscissa:  wt.«  PS.  Full  line,  open  symbols:  pjp.  Full  and  dashed  lines, 

full  circles:  Tg  p  and  p  (see  text).  Full  squares  (PS-d3)  and  tri¬ 
angles  (PS-dg):  Tg  SE-.  Dotted  line:  as  determined  by  dielectric  re¬ 

laxation  and  thermal  methods  (ref.  7). 
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^  A  theory  of  nuclear  magnetic  interaction  is  based  on  the 
study  of  the  stochastic  rotation  operator.  The  theory  is  applied 
explicitly  to  relaxation  by  anisotropic  chemical  shift  and  to  spin- 
rotational  interactions.  It  is  applicable  also  to  dipole-dipole  and 
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Theory  of  Nuclear  Magnetic  Relaxation 
James  McConnell 

Dublin  Institute  for  Advanced  Studies 
Dublin  4,  Ireland 


A  theory  of  nuclear  magnetic  relaxation  centred  around  the 
rotation  operator  is  proposed.  The  theory  is  applicable  in  principle 
to  molecules  of  arbitrary  shape,  and  account  is  taken  of  the  effects 
of  their  inertia.  Calculation  of  relaxation  times  associated  with 
anisotropic  chemical  shift  and  spin-rotational  interactions  provides 
illustrations  of  how  the  theory  may  be  employed. 
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1.  INTRODUCTION 

The  theory  of  nuclear  magnetic  relaxation  associated  with  random 

1  2 

thermal  motion  has  been  presented  by  Abragam  and  by  Hubbard  .  In 
the  present  investigation  a  theory  is  proposed  in  which  the  rotation 
operator  plays  a  central  and  explicit  part.  This  is  facilitated  by 
employing  results  for  the  rotational  Brownian  motion  of  a  rigid  body 
of  arbitrary  shape,  that  were  derived  by  calculations  in  which  the 

■3 

effects  on  the  inertia  of  the  body  are  included  . 

In  the  next  section  the  stochastic  rotation  operator  will  be 
defined,  and  an  outline  will  be  given  of  the  way  in  which  the  operator 
and  its  ensemble  average  are  calculated.  In  section  3  it  will  be  shown 
how  relaxation  times  can  be  deduced  from  spectral  densities  associated 
with  the  rotation  operator.  This  will  be  illustrated  in  section  4  by 
applications  to  relaxation  by  anisotropic  chemical  shift  and  by  spin- 
rotational  interactions. 


d&^ritA 


2.  THE  STOCHASTIC  ROTATION  OPERATOR 


s'  " gfj  The  concept  of  rotation  operator  may  be  introduced  by  taking 

^  a  set  of  rectangular  coordinate  axes  0%  t  Ofy@Z  and  a  function 

..  of  t2.)  of  the  coordinates  ^,^,2)  of  a  fixed  point  P  • 

c  Let  us  now  rotate  the  coordinate  axes  about  the  origin  to  Ox'  Oh!  Dz* 

D  *  a  t 

and  consider  the  same  function  of  the  new  coordinates  (x/,y./Jz') 

of  .  We  write 

r  'Rifh>v>2)  (1) 

and  we  say  that  is  the  rotation  operator  associated  with  the 

rotation  of  axes  and  with  the  function  ^  .  When  has  more 

than  one  component,  eq.  (1)  is  to  be  interpreted  as 


V- 


It  is  helpful  to  introduce  infinitesimal  generators  of 
rotation  Jx  ,  ,  J2  by 


Or  ^  we  rotate  Ox ,  Dy  fiz  through  an  angle  about 

.  .  .  ..  .  .  ^  ...  .  -  '  ,  1 


an  axis 


/.r  "-tt" 


0 

*■'■**■ 


*yr 


specified  by  the  unit  vector  *£  ,  it  is  easily  proved  that 

$  »  (  J‘  ^0-  (3) 

If  a  rigid  body  is  rotating  with  angular  velocity  U>[0)  about  ■€  , 
it  follows  from  (3)  that 


•f<fe 

*  -  V£/J!„  r°t 


The  rotation  of  the  coordinate  axes  may  be  achieved  by 


rotating  about  the  Z -axis  through  an  angle  V  ,  then  about  the 
-axis  through  an  angle  ji  and  lastly  about  the  ?  -axis 
through  an  angle  ^  **.  We  see  from  (3)  that  *|^  is  now  given  by 

R  -  <s> 


4 


••  .  If  for  this  rotation  we  apply  (2)  to  the  set  of  spherical  harmonics 

W  ***.  V  (  4).„e 

f  ..  .  X  ^<*-9  Cos  <j>  ,  y  -  ,  2  -  +00%$  y 

then  it  is  usual  to  write  (2)  as 

t 

V9/  ^ '  /* : /'j  Yjl  *  ft  ,6) 


T)>  J 

where  .  J/ .  ,  are  Wigner  functions.  Equation  (6)  shows  that 

m  j  9f\  n I 

J/  is  the  >?/tn  -element  of  the  matrix  representative  of  ^ 

in  the  representation  with  basis  elements  "Y* •  .V  , V* 

that  is  to  say  *  J  J>  J  JJ 


that  is  to  say> 


Xy>M  •  'R- 


We  apply  these  ideas  to  a  molecule  that  is  undergoing  steady 

Jt-c.apr  orr*tfi,  state  Brownian  motion.  Then  the  angular  velocity  lo(t)  in  (4) 

/  /Lj£f  i«rc  ) 

is  a  random  variable,  and  so  also  is  the  solution  of  this  equation. 


We  refer  to 


nit> 


as  the  stochastic  rotation  operator.  To 


obtain  information  about  CO 


we  take  axes  fixed  in  the  molecule 


with  origin  at  the  centre  of  mass  and  in  the  directions  of  the 
principal  axes  of  inertia.  We  assume  that  the  components  of  CJ 
referred  to  the  molecular  coordinate  system  obey  the  Euler-Langevin 
equations  ^ 


Up 
P.  <\  Cyf, 


'£££>  l-ft  "(I,- 1 

where  Xt  ,  Xx  t  Xj  are  tllc  principal  moments  of  inertia, 

£?.£>,  are  the  components  of  angular  velocity, T  4  7  A  7  A 

the  components  of  the  driving  couple  and 

*^the  components  of  the  frictional  couple 
resisting  the  motion.'  From  (8)  we  may  deduce  the  value  of  the  time- 
correlation  function  of  two  components  of  angular  velocity,  which 
will  be  useful  for  subsequent  calculations. 


«^e 

cp. 


In  order  to  solve  (4)  we  employ  a  method  of  solution  of 
nonlinear  differential  equaitons  that  goes  back  to  Krylov  and 
Bogoliubov  ^ .  It  is  assumed  that  the  solution  Rio  consists  of 
a  slowly  varying  ensemble  average  >  about  which  there  are 

random  fluctuations,  and  the  solution  is  expressed  as 

( 1  ♦  t  F (t)  +  ^  FTo  +  •  •  ■ )  <  R  to)>,  (9> 

where  J  is  the  identity  operator  and  i  is  a  small  dimen¬ 
sionless  parameter.  ■R(t)  describes  the  rotation  of  the  molecule 
from  its  orientation  at  time  zero.  The  ensemble  average 
is  nonstochastic  and  it  is  supposed  to  satisfy  some  equation 

~~  o) 

By  using  the  knowledge  that  we  have  obtained  from  (8)  we  can  often 

deduce  from  (4),  (9)  and  (10)  the  values  of  i  /t) 

r\(n  tow  *  1  i 

'  for  substitution  into  (9)  and  (10).  We 

may  then  be  able  to  solve  (10)  for  •  Analytical 

solutions  have  in  fact  been  found  for  a  molecule  that  is  spherical 

-} 

or  linear  or  a  symmetric  top  or  even  an  asymmetric  top  .  For 
most  NMR  problems  the  value  of  is  all  that  we  need.  If 

we  do  need  the  value  of  Rid)  ,  we  can  obtain  it  from  (9). 
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*  3.  BASIC  THEORY  OF  NUCLEAR  MAGNETIC  RELAXATION 

| 

To  study  nuclear  magnetic  relaxation  we  take  a  set  of 
cartesian  axes  fixed  with  respect  to  the  laboratory.  A  constant 
(.{ytJ  strong  magnetic  field  H0  acts  in  the  2  -direction.  If  f-fg 

„  -  interacts  with  a  nucleus  of  spin  J  and  gyromagnetic  ratio  Y  , 

Z  _n_  1 1  7  /  Qt> 

the  field  produces  a  Hamiltonian  —  y IfZ  J2  ,  denoted  also  r 

*y  *  &>  >  with  non-degenerate  energy  levels.  We  now  use  the 

*  normalized  eigenfunctions  of  ax  as  a  basis  to  express  matrix 
elements  of  operators.  r^^»ffi^perturbing  Hamiltonian  a  G(t) 

A 

like  that  arising  from  dipole-dipole,  quadrupole,  chemical  shift 
or  spin-rotational  interactions  causes  a  relaxation  process,  and 

*  the  results  are  usually  expressed  in  terms  of  longitudinal  and 
transverse  relaxation  times  "Tj  and  T,  ■  T  referring  to  the 

2  -direction.  For  any  particular  interaction  one  has  two  problems: 
a.  to  express  ”T,  and  in  terms  of  spectral  densities, 

'  b.  to  calculate  the  spectral  densities. 

.  The  most  direct  way  of  handling  a.  is  to  employ  a  general 

X.  r.  0C  Q  KA 

method  due  to  Redfield  .  If  with  -T*  "K,  y ,  Z  is  a 

component  of  magnetic  moment  per  unit  volume  due  to  the  inter- 
i  t.  V  V°  action  and  PAA*  denotes  an  element  of  the  density  matrix, 


-t'f- 


method  due  to  Redfield 


c.  «yr  "-ro  action  and 
then 

■t.r.  op  'CttX 

J?-  v -tft . 

where 


d.<MJ 

dt 


S'  V.  -  rJ  fol' ~  J  (W, -W,) 

6,  J*n  rN-w'A»'  *  r 


y.r.tp 
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The  J  (<*}  )  is  a  spectral  density  referred  to  in  a.  and  b. 

C  is  an  abbreviation  of  (  (A  I  >  where  \<A 

are  eigenfunctions  of  £  tyo  .  The  energy  of  the  state  \  &) 
is  £.  tPj,  .  On  substituting  (12)  and  (13)  into  (11)  we  may 
obtain  an  equation 


dOV> 

ctt 


From  this  would  follow 
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To  see  how  the  spectral  densities  are  to  be  calculated  ^  ^ 
we  shall  first  of  all  take  a  special  form  of  G(t),  namely, 


Z  -  \ 


where  Fjt)  #  F,  (t)}  ?  fjt)  , 

rotations  like  the  spherical  harmonics 


transform  under 


J.f _  y 


,  IP't  ' 

Cm 


*-)*  »  ■  p- 

This  will  permit  the  application  of  (6)  and  (7)  to  Fit)  •  We 

^  f 

have  from  (16) 

^  •  "i  /-)*'  f^w  1  a, 

£r-t  ' 

where  is  the  adjoint  of  the  operator  deduce  that 
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The  angular  brackets  denote  here  an  average  over  the  random 
variables,  that  is  to  say,  the  angular  velocity  components, together 
with  an  average  over  the  initial  orientations  of  the  molecule.  A 
study  of  the  transformation  properties  of  spherical  harmonics  shows 

Q 

that 


'V 


/H' 


-*,>/:  -2 

a  constant,  and  <  R(w>  is  the  ensemble 


+*>€nU,  *.<L 


where  F/^F^ 
average  of  the  rotation  operator  for  the  molecule  in  question.  We 
see  that  we  may  write 


independent  of  ^ 


where 


,  and,  from  (13)  and  (17),  that 


1* .  *  i  z"  f; 

_  -O*  1 

•’i.w  * 


(18) 


(19) 


We  know  K 
V 


from  the  nature  of  the  perturbing  Hamiltonian 
that  gives  rise  to  the  interaction  and  we  know  ^  JP^^^from  the 
investigation  of  the  rotational  Brownian  motion  of  the  molecule. 
Hence  we  can  in  principle  calculate  JVw)  from  (19)  and 
consequently  the  relaxation  times  from  (11),  (12), (13), (15)  and 
(18).  Moreover,  since  the  integral  in  (19)  is  independent  of  the 
interaction  Hamiltonian  t.  CM  ,  results  for  relaxation  times 
can  easily  be  taken  over  from  one  interaction  satisfying  (16)  to 
another  by  replacing  the  values  of  in  (19)» 
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4.  CALCULATION  OF  RELAXATION  TIMES 
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As  applications  of  the  Redfield  theory  we  shall 
investigate  two  types  of  interaction  that  cause  nuclear  magnetic 
relaxation.  The  first  is  anisotropic  chemical  shift,  which  is 
an  example  of  an  interaction  satisfying  (16).  The  other  is  spin- 
rotational  interaction,  which  does  not  obey  (16). 

For  anisotropic  chemical  shift  it  may  be  shown  that10, 

A0-  mi*  ,a,=  -4iu,a,- 4bRi 


where 


and  that 


r-i/i  . f; . o .  K-ltK-,  • 

where  the  non-vanishing  elements  of  the  diagonalized  anisotropic 


'  K- 

S' 

Pi  \U 


chemical  shift  tensor  are 


"a  (  J)  4'  >  1  1  P  ,  K‘  ' 


Hence,  from  (l8), 


'  >) r  i  J„  <<*  { *  I  At  lp)(  /S' I  A*  I  °< ' ) 

and  Jc$(^  is  given  by  (19)  with  the  values  of  F  taken  from 
(21).  Since 

<M>r*/v<u 

where  (V  is  the  number  per  unit  volume  of  nuclei  taking  part  in  the 
relaxation  process,  we  may  obtain  the  relaxation  times  by  applying 
(11),  (14)  and  (15)  to  .  Employing  (20)  we  find  after  an 

elementary  calculation  that  10 


10 
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'S  T  3  Hp  Jcs^o) 


is  =  [  h  i5  f°J  +  ^  <J ffu]. 

»  7- 


(22) 


To  complete  the  calculations  we  must  obtain  an  express¬ 
ion  for  fw)  from  (19)  and  (21).  What  was  equivalently  the 

same  thing  was  done  by  Abragam  but  only  for  a  spherical  molecule 

12 

and  with  inertial  effects  neglected  .  The  explicit  use  of  the 
rotation  operator  has  allowed  us  to  include  the  inertial  effects 
and  to  extend  the  discussion  to  a  linear  molecule  and  to  a 
symmetric  top  molecule  7  We  shall  put  down  the  result  only  for 
the  latter  case  taking  the  nucleus  to  be  at  the  centre  and  the  axis 
of  symmetry  as  the  third  coordinate  axis  in  (8),  so  that  J,  -  J 

3  ' 1 


j£Sm-  5^r  t  (n.+AD,) 
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where 


4T 
L  h 


(23) 


The  relaxation  times  follow  from  (22). 

We  see  from  (16),  (18)  and  (19)  that,  when  the  interact¬ 
ion  Hamiltonian  can  be  expanded  as  a  linear  combination  of  elements 
of  a  five-dimensional  spherical  tensor,  we  require  only  ^^/tJ^for 
the  molecule  in  order  to  calculate  the  spectral  density  J(u>)  and 
consequently  the  relaxation  times.  This  procedure  fails  in  the 
study  of  nuclear  magnetic  relaxation  by  spin-rotational  interactions. 


* 'ft 


Then 


tC-t(l-Cj), 


(24) 


where  1  is  the  spin  operator  of  the  nucleus  that  interests  us, 
is  the  three-dimensional  spin-rotational  tensor  and  ^  J  is 
the  angular  momentum  operator  of  the  molecule  that  contains  the 
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jay 


nucleus.  Since  the  number  of  independent  elements  of  C  is  not 

five,  the  of  (24)  cannot  be  expanded  as  in  (16). 

A  general  theory  of  nuclear  magnetic  relaxation  by  spin- 

rotational  interactions  has  been  based  on  the  stochastic  rotation 

operator  The  angular  momentum  J  is  replaced  by  its 

classical  value  (~l  u  T  w  7  tu. 1  ,  We  have  found  that  the 
v  x,  ^  sr  s-r 

Redfield  method  yields  relaxation  times  I  /  given  by 

•  >  la 


m ,  or  —  s-r  =  ^  } 

where 

/)  vr 


where 


./•<•  V  "t*** 

*  r  "1*K 


jJ  <27) 


/U/V-f 


-  6^  ^  t  6 


and  denotes  matrix  representation  with  respect  to  the  basis 

•>*  10  li  '  ^  v 

We  see  from  (25),  (26),  (27)  that  a  knowledge  of  \  JvfrV 
is  no  longer  adequate  for  the  calculation  of  spectral  densities  and 
relaxation  times.  We  must  now  use  our  knowledge  of  the  values  of 
J-  r.  <y 1.  tF7w,t  and  of  to  calculate  'R/tj  from 

(9).  Then  we  must  obtain  ^  ^ It)  ^  IcJuj^Io/ >  calculate 
its  Laplace  transform  and  the  matrix  elements  of  this,  evaluate  CfS ) 
from  (27)  and  employ  (25)  and  (26)  to  deduce  the  relaxation  times. 

This  has  been  done  in  principle  for  an  asymmetric  molecule  and  with 
the  inclusion  of  inertial  effects.  Explicit  results  can  be  given 
for  molecules  that  are  linear,  axially  symmetric  or  spherical. 

For  simplicity  we  shall  report  these  results  only  in 
the  extreme  narrowing  approximation  of  £uc  <i<.  ~T / (lj})  <~mp 

-r  9*  «  .  4^  ’’ 

where  I  and  r?  stand  generically  for  the  moments  of  inertia  ,  .  „ 

-r-  S-r  -r  St  J 

and  the  friction  constants  in  (8).  Then  and  f  are  equal, 

and  we  shall  denote  their  common  value  by  ""H,  •  In  the  case  of 
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1  1 

a  spherical  molecule 

J--  2W  [//•'. 


f.t  Cf,  "'blfa*  where  K  -  ^£T /{[%*)  ^the  spin-rotational  tensor  component  Cn 

((  ^ 

.  refers  to  the  direction  of  the  radius  through  the  nucleus  and 

refers  to  a  perpendicular  direction.  In  the  case  of  a  molecule 
which  has  a  principal  axis  of  inertia  through  the  centre  of  mass 
that  passes  through  the  nucleus  and  which  is  such  that  this 
i  principal  axis  is  an  axis  of  symmetry  (_  for  the  molecule 

With  X*  ,  3 

-±r  ,  ill  (  -  ,  2,  C 

T„  +  B,  +  ST>, 


+  *T|V 5’ +  fvi,' 

In  this  equation  ^/|  >  ^3  >  ?5  refer  to  the  axis  of  symmetry, 

C ,  i.  >  3,  refer  to  a  perpendicular  axis  and  D,  >  T^, 

are  defined  in  (23).  In  the  case  of  a  linear  molecule 


,  /„2  (29) 

c\ 
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Calculations  have  been  performed  for  the  above 
molecular  models,  when  inertial  effects  are  ignored  It  has 

been  found  that  (28)  is  altered  to 

_l  nrr(C +2C) 

T„  '  It'S  ""  ’ 


that  (29)  is  altered  to 


J_  z£T 

T„  31' 


n,C  .  iC 

'T'v  '  _ 


and  that  (30)  is  unaltered.  The  inclusion  of  inertial  effects 


produces  for  the  sphere  and  for  the  symmetric  rotator  corrections 
of  order  K  ,  which  is  at  most  a  few  per  cent.  There  is  zero 
correction  for  the  linear  molecule,  even  though  the  calculations 
have  been  performed  to  an  accuracy  of  K*  . 


5-  CONCLUSION 

An  analytical  theory  of  nuclear  magnetic  relaxation  has 
been  based  on  the  stochastic  rotation  operator  and  on  the  Redfield 
method  of  establishing  a  differential  equation  for  the  ensemble 
average  of  the  magnetic  moment  produced  by  a  specified  interaction 
mechanism.  The  theory  has  been  applied  to  examples  of  two  different 
classes  of  interaction,  namely,  anisotropic  chemical  shift  and  spin- 
rotational.  The  present  theory  is  confined  to  the  study  of  the 
behaviour  of  a  single  molecule. 
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From  the  very  earliest  days  of  NMR  it  has  been  possible  to  study  the 
motion  of  molecules  in  the  solid  state,  first  from  the  narrowing  of  the  proton 
dipolar-broadened  spectrum,  and  in  more  detail  from  the  measurements  of  the 
proton  spin-lattice  relaxation  time  Tj.  So  we  wondered  how  much  we  could 
learn  about  the  dynamical  behavior  of  larger  molecules  such  as  proteins  from 

* 

NMR  in  the  solid  state.  This  has  the  interest  that  we  are  examining  the 
molecules  in  the  same  environment  as  that  in  which  their  X-ray  structures  were 
determined  and  it  would  also  complement  NMR  studies  on  these  molecules  in 

solution. 

i 

One  advantage  of  working  in  the  solid  state  is  the  very  wide  temperature 
range  which  is  accessible  from  room  temperature  down  to  quite  low  temperatures 
and  this  gives  the  chance  to  characterize  the  motions  well  with  activation 

j  parameters.  In  solution  often  only  a  small  temperature  range  is  accessible 

before  denaturing  and  degradation  occur. 

A  second  advantage  of  working  in  the  solid  state  is  that  the  large 
protein  molecule  is  essentially  fixed  in  the  crystal  lattice  so  that  we 

» 

»  examine  directly  the  intramolecular  motions.  In  solution  on  the  other  hand 

the  intramolecular  motions  are  superposed  on  tumbling  of  the  whole  molecule, 
which  complicates  the  analysis,  and  of  course  it  can  only  be  applied  to  those 
groups  whose  spectral  lines  are  separately  resolved  and  assigned. 

4 

In  the  present  work  we  are  looking  at  the  relaxation  of  the  whole 
collective  assembly  of  dipolar-coupled  protons.  So  we  cannot  be  specific 
about  the  details  of  particular  individual  groups,  but  as  we  shall  see  we  can 
determine  the  types  of  motion  which  contribute  most  and  we  can  characterize 

i 

them.  A  protein  contains  several  hundred  protons,  maybe  a  thousand  or  more, 
each  with  its  own  molecular  environment,  which  makes  it  a  complicated  system 
to  investigate.  So  the  philosophy  of  our  work  was  to  start  with  simpler 
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precursors  and  having  established  some  types  of  behavior  and  understood  them 
quantitatively  we  could  then  move  on  to  protein  molecules  themselves  with  more 
confidence. 

So  we  began  by  looking  at  the  proton  relaxation  time  Tj  of  all  the 
classical  20  amino  acids  of  which  proteins  are  composed,  in  the  solid  state 
from  100  to  500K.  This  work  has  been  reported  In  a  series  of  three  papers 
(1,2,3).  Analysis  using  the  well-known  Kubo-Tomita  theory  of  dipolar 
relaxation  (4)  developed  from  the  original  relaxation  theory  of  Bloembergen, 
Purcell  and  Pound  (5),  quantitatively  identified  reorientation  of  methyl  and 
NH^  groups  and  motions  of  other  side  groups.  At  each  temperature  a  single 
correlation  time  characterized  each  motion  and  its  temperature-variation 

followed  a  simple  Arrhenius  activation  law,  from  which  an  activation  energy 
and  pre-exponential  factor  was  determined.  The  whole  relaxation- 

temperature  curve  was  fitted  by  computer,  not  just  the  asymptotes. 

Measurements  were  made  at  two  frequencies  and  the  same  parameters  were  used  in 
fitting  the  results  from  both.  A  comparison  of  the  measured  and  calculated 
dipolar  coupling  constants  provided  a  strong  test  of  the  correct 
identification  of  the  motions.  This  work  on  the  solid  amino  acids  is  a 
subject  in  itself  and  some  of  the  amino  acids  exhibited  most  interesting  side 
chain  motions,  particularly  arginine,  cysteine,  methionine,  phenylalanine  and 
valine. 

These  measurements  were  extended  to  polycrystalline  dlpeptldes  and 
tripeptides  (6),  for  example  alanylglycine  and  glycylalanlne  and  the  series 
glycine,  diglycine,  triglycine,  polyglycine.  We  usually  found  the  features  of 
the  monomeric  amino  acids  carried  over  with  minor  modifications  into  the 
peptides.  Next  we  examined  a  series  of  homopolypeptides  (7)  and  by  then  felt 
we  had  enough  experience  to  help  us  tackle  some  solid  proteins. 

The  four  proteins  examined  in  polycrystalline  form  were  a-chymo trypsin, 
insulin,  lysozyme  and  ribonuclease  A,  selected  because  they  are  well- 
characterized,  their  X-ray  structures  have  been  determined,  they  are  readily 
available  in  pure  form  and  are  not  too  expensive.  We  used  about  0.5g  for  each 
specimen,  and  pumped  them  for  24  hours  at  room  temperature  to  remove  oxygen 
and  most  of  the  water.  Measurements  were  made  at  18,  30  and  60  MHz  on  each 
protein  from  300K  down  to  10K.  Results  for  ribonuclease  A  are  shown  in  Fig.  1. 

The  results  (8-11)  exhibited  altogether  broader  shallower  curves  than 


Chose  found  for  Che  monomeric  amino  acids  and  Che  simple  pepCides,  and  could 
only  be  accounCed  for  assuming  a  disCribucion  of  correlaCion  cimes  This 

is  really  noc  surprising.  In  Che  pure  amino  acids  we  somecimes  resolved  cwo 
or  even  chree  correlaCion  cimes  from  cwo  or  Chree  independenc  molecular  group 
mocions  ac  each  Cemperacure.  In  a  biopolymer  such  as  Chese  protein  molecules, 
with  very  many  different  proton  sites  and  degrees  of  freedom,  many  correlation 
times  are  needed.  We  should  note  that  we  could  fit  any  one  of  Che  relaxation 
curves  measured  at  a  single  measuring  frequency  with  a  single  thermally- 
activated  correlation  time,  but  not  all  three  curves  for  the  three  measuring 
frequencies  simultaneously.  This  illustrates  how  important  it  is  to  work  at 
several  measuring  frequencies. 

Four  well-known  distribution  functions  were  tried:  Gaussian  (log¬ 
normal),  Cole-Cole,  Cole-Davidson,  Fuo as -Kirkwood.  It  was  found  tht  the 
Gaussian  distribution  provided  the  best  characterization  of  the  data  for  all 
four  polycrystalline  proteins  over  the  range  70-250K. 

We  can  go  a  long  way  towards  identifying  the  motions  responsible.  The 
earlier  work  on  precursors  leads  us  to  expect  methyl  groups  and  their 
reorientation  to  be  important  in  the  side  chains  of  six  amino  acids,  namely 
alanine,  isoleucine,  leucine,  methionine,  threonine  and  valine.  The 
relaxation  minima  of  the  proteins  occur  in  approximately  the  same  temperature 
range  as  in  most  precursors.  The  activation  energies  are  very  similar  to 
those  of  the  amino  acids;  we  do  not  expect  them  to  be  exactly  the  same,  but  it 
is  satisfactory  that  they  are  similar  and  it  suggests  that  most  of  the 
restrictions  to  reorientation  are  intramolecular.  Moreover  we  can  calculate 
the  expected  relaxation  constant,  assuming  that  rapid  spin  diffusion  maintains 
a  common  spin  temperature  and  that  the  methyl  proteins  have  to  relax  all 
protons  in  the  protein  molecules.  This  analysis  shows  that  methyl 
reorientation  In  the  side  chains  accounts  for  about  702  of  the  measured 
relaxation  constants,  leaving  some  302  to  be  accounted  for  by  side  chain 
motions,  segmental  motion,  reptation  and  whole-body  motions. 

Since  we  are  dealing  with  an  inhomogeneous  system  with  methyl  and  other 
groups  in  a  wide  variety  of  environments,  the  assumption  of  a  distribution  of 
values  of  xc>  each  characterizing  an  exponential  correlation  function,  seems 
physically  more  appropriate  than  use  of  a  single  suitably-chosen  non- 
exponential  correlation  function  such  as  that  of  Williams  and  Watts  (12). 


Above  250K  and  below  60K  the  experimental  values  of  Tj  fall  below  the 
calculated  curves  for  all  four  proteins  and  evidently  other  processes  become 
important  here.  At  low  temperatures  there  seems  to  be  a  weaker  but  less- 
hindered  process.  To  cast  light  on  this  we  examined  two  homopolypeptides 
devoid  of  methyl  groups,  namely  polyproline  and  polyglycine. 

Polyprollne  showed  a  surprisingly  strong  proton  relaxation  (13). 

Analysis  of  the  relaxation  constant  showed  that  this  could  be  attributed  to 
conformational  motion  of  the  proline  rings.  There  are  a  number  of  prollne 
residues  in  proteins,  for  example  there  are  nine  in  a-chymotrypsin.  This 
motion  can  therefore  be  expected  to  contribute  significantly  to  proton 
relaxation  over  the  whole  temperature  range. 

By  contrast  the  proton  relaxation  in  polyglycine  is  much  weaker  (13), 
arising  from  segmental  motions  of  the  chains  which  modulate  the  dipolar 
interactions  of  the  methylene  and  imlno  groups.  Below  60K  the  relaxation  of 
solid  proteins  and  polyglycine  tend  to  become  equal;  this  supports  the  view 
that  segmental  motions  of  the  protein  chain  make  important  contributions  to 
relaxation  in  solid  proteins  at  low  temperatures. 

Above  250K  the  deviation  of  Tj  below  the  computed  curves  is  attributed  to 
additional  relaxation  contributions  of  water  molecules  in  the  structure.  We 
have  investigated  this  with  measurements  on  polycrystalline  lysozyme  and 
a-chymotrysin  progressively  hydrated  with  H20  and  D20  up  to  60%  by  weight 
(14).  The  amount  of  water  added  was  measured  by  increase  of  weight  and  by 
Karl-Fischer  titration.  The  Karl-Flscher  titration  showed  that  our  normal 

specimens  pumped  at  room  temperature  retained  about  2 %  water,  strongly  bound. 

S«  foj.X. 

Two  main  features  were  noticed  (14),A  Above  180K  increase  of  hydration 
progressively  decreased  Tj  exhibiting  a  characteristic  minimum  which  shifted 
to  progressively  lower  temperatures.  Below  180K  Increase  of  hydration 
progressively  Increased  Tj.  The  increased  relaxation  rate  above  180K  is  due 
to  reorientation  of  the  water  molecules.  Additional  water  is  bound  less 
strongly  causing  the  Tj  minimum  to  shift  to  lower  temperatures.  Analysis 
shows  (14)  that  the  activation  energy  for  this  motion  decreases  from  30kJ/mole 
for  the  first  increments  of  water  to  20kJ/mole  for  the  last  Increments.  It  is 
noteworthy  that  the  last  increments  of  water  contribute  to  relaxation  even 
down  to  180K.  The  relaxation  constant  for  the  water  reorientation  turns  out 
to  be  about  56%  of  that  for  isotropic  intramolecular  dipolar  relaxation  of 


pure  water.  This  reduced  value  can  be  attributed  to  three  causes:  (1) 
anisotropy  of  the  reorientation  of  the  water  molecules  on  their  sites,  (2)  a 
distribution  of  activation  parameters  on  the  various  sites  of  attachment,  (3) 
a  lengthening  of  the  interproton  distance  due  to  hydrogen  bonding  of  the  water 
molecules  on  their  sites. 

Below  160K  Tj  increased  uonotonlcally  and  approximately  linearly  with 
successive  additions  of  water.  At  these  lower  temperatures  the  water 
molecules  are  frozen  in  the  structure  and  contribute  nothing  to  relaxation, 
but  they  do  add  to  the  load  of  protons  needing  to  be  relaxed  by  other  protein 
motions. 

When  D2O  is  added  instead  of  HjO  we  have  similar  effects,  but  very  much 
weaker  because  of  the  relative  inefficiency  of  deuterons  in  relaxing  the  load 
of  protons;  indeed  calculations  shows  that  it  is  U.2X  smaller  which 
approximately  accounts  for  what  is  observed. 

This  study  confirms  that  the  presence  of  water  in  proteins  causes  extra 
relaxation  at  higher  temperatures  due  to  motions  of  the  water  molecules 
themeselves,  but  there  16  no  evidence  of  any  change  In  the  protein  molecules 
induced  by  the  presence  of  the  water. 

All  the  work  described  so  far  Is  based  on  measurements  of  T^  which  are 
particularly  responsive  to  dynamic  behavior  in  the  frequency  range  105  -  1010 
Hz.  Recently  we  have  extended  our  studies  to  measurements  of  the  dipolar 
relaxation  time  Tjq  using  the  Jeneer-Broekaert  pulse  sequence  (15),  since  this 
is  responsive  to  dynamical  behavior  in  the  spectral  region  1  to  10^  Hz, 
enabling  a  wide  spectral  range  1  to  10*®  Hz  to  be  investigated  overall. 

Extending  the  expression  for  T^  for  a  distribution  of  correlation  times 
in  the  same  manner  as  we  did  previously  for  T^  and  using  the  activation 
parameters  previously  obtained  from  measurements  of  T^ ,  we  predict  that  a 
minimum  in  Tjp  should  be  found  at  about  85K,  about  100K  lower  than  for  Tj, 
Experiment  bore  this  prediction  out  closely  for  both  a-chymotrypsln  and 
lysozyme  (16).  This  extrapolation  was  a  severe  test  of  the  parameters 
obtained  from  Tj  measurements  and  gave  some  confidence  to  the  description  and 
characterization  of  the  motions  obtained  through  Tj.  The  shift  of  the  TjD 
minimum  100K  lower,  where  the  motions  are  some  10*  times  slower  gives  the 
opportunity  to  examine  new  motions  at  higher  temperatures.  The  minimum  value 


of  TjD  is  210  times  shorter  than  that  of  Tj  at  60MHz  for  lysozyme,  but  theory 
predicts  a  ratio  of  3000.  A  similar  discrepancy  was  found  for 
a-chymotrypsin,  and  is  attributed  to  a  spin-diffusion  bottleneck. 

In  some  cases  the  use  of  enables  relaxation  mechanisms  and  motions  to 
be  resolved  which  are  not  resolved  in  Tp  An  example  was  provided  by  solid 
polyvaline.  In  addition  to  a  deep  minimum  in  attributed  to  methyl  group 
reorientation  in  the  monomer  side  chains,  there  is  a  weaker  minimum  at  about 
200K  which  is  attributed  to  a  more  hindered  reorientation  of  the  whole  side 
chain  which  modulates  a  much  smaller  fraction  of  the  overall  dipolar  second 
moment .  v -- 

x> 

So  summarizing,  At  has  been  possible  through  the  study  of  proton 
relaxation  in  solid  proteins  and  related  molecules  to  identify  and 
characterize  the  following  molecular  motions  which  contribute  to  relaxation: 

-Ml)  methyl  group  reorientation; 

^(2)  segmental  motion; 

—*■(3)  side  chain  motion  , 

»(4)  proline  ring  puckering j  a  ,  f 
(5)  water  molecule  reorientation.  . 

Finally  it  is  a  pleasure  to  express  my  great  Indebtedness  to  my 
colleagues  at  Nottingham  without  whom  these  studies  could  not  have  been 
undertaken,  especially  D.  J.  Bryant,  E.  M.  Cashell,  Q.  A.  Meng,  R.  Gaspar,  D. 
N.  Bone  and  T.  Z.  Rizvi. 
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Figure  Captions 


Fig.  1.  The  temperature  dependence  of  the  proton  spin-lattice  relaxation  time 
Tj  in  solid  ribonuclease  A  between  10  and  300K.  The  full  lines  are 
theoretical  curves  calculated  as  described  in  the  text. 

Fig.  2.  The  variation  with  temperature  of  the  proton  spin-lattice  relaxation 
time  T^  at  60  Mlz  is  polycrystalline  lysozyme  with  various  degrees  of 
hydration.  Full  circles  and  full  lines  refer  to  hydration  with  f^O; 
open  circles  and  dashed  lines  refer  to  hydration  with  I^O. 
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